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ABSTRACT 
Pewee Ib 


Methods have been developed for the determination 
OL MECrOSCOpICcrana Macroscopic acid dissociation constants 
of polyprotic acids from nuclear magnetic resonance 
chemical shift measurements. Depending upon the position 
of a nmr-active nucleus relative to two simultaneously 
deprotonating functional groups, its chemical shift may 
rerlect ithe stace of protonation of one Of the two groups 
(a unique resonance) or the state of both groups (a 
common resonance). The change in chemical shift of a 
unique resonance with changes in pH gives directly the 
fractional deprotonation of one of the two simultaneously 
deproconating Groups as 7a function of pH. wim order to 
calculate fractional deprotonation values from common 
resonance data, it is necessary to determine the total 
change in the observed chemical shift due to the depro- 
tonation of each of the two groups. Nonlinear least 
squares curvefitting techniques have been developed to 
Calculate both macroscopic and microscopic acid dissoci-— 
ation constants from fractional deprotonation data. Also, 
a nonlinear technique has been developed to calculate 
macroscopic constants directly from common resonance chem- 
ECavesh@it data, 

Using these methods, microscopic and macroscopic 


acta dissociation constants were calculated from proton 
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magnetic resonance chemical shift titration curves for 
the amino acids lysine and ornithine, the dipeptides 
lysylglycine and glycyllysine and the related compound 
ethylenediaminemonoacetic acid. Both macroscopic and 
Microscopic constants were calculated from carbon-13 nmr 
Litheat von Ounves (for 2,3-dlaminoproplonic cid) ornithine 
and lysine. The carbon-13 data for 2)/4-diaminobutyric 
acid proved to be anomalous and provided only macroscopic 
Gonseanis. 

In connection with the carbon-13 nmr studies of 
the amino acids, linear parameterization schemes have been 
developed to predict the cmr chemical shifts of the alkyl 
and carboxyl carbons of amines, carboxylic acids, amino 
acids, and related compounds in aqueous solution. The 
parameters were determined by multiple, linear regression 
analysis of the observed chemical shifts of forty-five 
model compounds. The chemical shift parameters should 
be of value in the assignment of observed chemical shifts 
to the carbons of peptides and proteins as well as amines, 


carboxylic acids, and amino acids. 
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The aqueous solution chemistry of trimethyllead 
perchlorate and of the trimethyllead complexes of five 
Carooxylic acids of PK, Values ranging from 3740) 16.47.95 
has been studied by proton magnetic resonance spectroscopy 
and potentiometry. Pmr is shown to be a direct method 
for studying the solution chemistry of the trimethyllead 
cation; the chemical shift of methyl protons and, to a 
lesser extent, the coupling constant for spin-spin 
coupling between lead-207 (I = 4%, natural abundance 21.1%) 
and the protons of the methyl groups are sensitive to the 
nature of the ligands bonded to the trimethyllead and 
thus provide information at the molecular level. The 
equilibrium constant. for the reaction of (CH,) ;Pb™ with 
hydroxide to form (CH,) ,PbOH has been determined from the 
pH dependence of the chemical shift of the methyl protons. 
Formation constants for the carboxylic acid complexes also 
were determined from the pH dependence of the chemical 
shift of the methyl protons in solutions containing both 
trimethyllead perchlorate and carboxylic acid. The form- 
ation constants of the trimethyllead complexes are found 


to anerease linearly as the acid dissociation constants 


Or the carboxylic acids decrease. 
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CHAPTER I 
INTRODUCTION 


The jactdabecse Chemistry of polyprotiac molecules 


is 


generally described in terms of stepwise, macroscopic acid 


dissociation constants. -.For-example, the titration of 
Lully protonated form of the amino acid glycine yields 
HNCH.CO.H 
5 Dien 2 
Macroscopic acidity Constants, pK) = 2.3 and pK, = 9.8 
(1). The macroscopic deprotonation scheme can be des- 
cribed by: 
+ + 
HOA Pee aA (1) 
+ 
K, = [H ey (2) 
[HA] 
HA + H + A™ (3) 
<— 
See Hee | A 
Ri aac in tas a: 


On the molecular level, the monoprotonated form HA can 
be represented by two microscopic forms as iidustrared 


Scheme I. 
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BY analogy tO acetic acid (pK, = 4,76 (2)) and ethylamine 
(pK, = 10.61) (3 )))e- pK, of glycine can be assigned to the 
-CO5H group and pk, to the amino group. Due to the widely 
Varying acadities of the two functional groups, the de- 
protonation will occur almost entirely via the upper path- 
Way eOumoCHeme Meanhdyveryeluttle OL the olycaine waillvbe 
DEeseVve as PCO ltl ateeany. pil. 

When the strengths of the acid groups are similar, 
two or more functional groups may simultaneously deproton- 
ate over the same pH range, and the macroscopic constants 
GObtaiticd by puatiltracionscan no longer be assigned Go a 


specific group. The amino acid cysteine, 


HSCH,CH-CO5H 


+ 
NH, 


which has been studied quite intensively, is a molecule 


Of tiis type. whe macroscopic, deprotonation ‘scheme of 
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cysteine can be represented as follows 


NS geen ep gee aha (5) 
3 <—— 2 
- fs 
HA —— HA + #H (6) 
rei 
= = + 
HA —— “-~ A +H (7) 
eerhee 
where PR) = Lei ay PK, = S33 and pK, = LO wi Be 44) 7. PK} 


Can’ be assigned to the carboxylic acid group but the 
assignment of pK, or pK to either the sulfhydryl or amino 
group, by analogy to similar Senses, is ambiguous. 
tTAtecially, Conn and Edsall (5) assigned the higher pk to 
the sulfhydryl group, but this assignment was reversed by 
Galvin (6)0. cudsel ly, as citedeby Ryklanvandeschmide ((7/)), 
essumed sthe intrinsic proton abiinities of the sulthydry! 
and amino groups, to be Similar, and thus considered the 
macroscopic pK values to be due to the deprotonation of 
DOLHECrOUDS  winic latter interpretation has proven to 

be correct and the microscopic acid dissociation scheme 


for the second and third deprotonations of cysteine is 


PAlGstravcea ain Scheme Ll (8). 


din AE = Gta ne 8 = oy a 


en? sa Giese fice “>. flames > WY ae cotuee og 
onlite 3 Key divs égg era y hs ot NG 31h. sad 7a ‘innate 
“* ay ¢ — 4 


ae 
aeeupidan 4h  ahnuqyhirs Phwig ©+. “pT ad&. Wet ‘SeOTe 
| a 2 + _ 

os Sf sstp kt nat Gonginae ‘c lp ahs’ One saw alta tate» 
(“leap 1 gay yacht, 2 


ar? a. ot ot! 4 rh 


ya Pe yee e maeyyy a} 7. hee | 
AT) abiedie’ bed nehiyn vo bento eh) aliditt. 4 1 a aie 
ht ‘ 


sa i : Ree ; 
Tesi ine os 6 AS04001 14 coset trentiitgear say ue 


eds (tareheinhed dudy Gas, c6d isa ad os BOI mr bof ae 
; ’ 2 he - 
ae aT | oe tp nee 
ets) Od, Tar Myregtie a Bj aa ngef ae Oe Fi : aRsere Pe ond Phy S20ns3a 

| ‘ a. Z 


OF, AHVSS GAY LRAT OI. Pasa. Tak eit 


engine tons Hb olege:) 6 bisa sim: clieadl af site sees, 


oer: 
Se 
v4 H,C-S Sse 
(tay) 
al ie 
(Gy) Pome eee (IV) 
H.,C-SH HC-S™ 
(2 
HC-NH., 
H.C-SH 
(nt) 
Scheme ITI 


The microscopic constants are denoted by Vower case ks 


and are defined by Equations 8 through ll. 
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The acid dissociation step to which a given microconstant 
refers is indicated by the subscript, the last number de- 
noting the group undergoing deprotonation and the preceding 
number(s) indicating groups which have previously under- 
gone deprotonation. The carboxyl, sulfhydryl, and amino 
GEOUpssare, denoted by the subscripts (1) (2), and (3) 
respectively. The macroscopic constants can be defined 


on the molecular level by: 


oF - + 
sae altet UVES ed ns ell sec Giese ae gly ETI) 
oo er Tae Ooi, eee vey 
= = + 
ee SESS tee Details bel 
K = (is) 
= ae ea ae 
5} [HA] 1 Pa came 


By subsettution of Rquations 8—ll into Equations 12 and 
sy eecan be shown that. the macroscopic constants are 


composites of the microscopic constants. 


K eis sie + k (14) 
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Top tee 
From the: preceding equations; it can ‘be seen that in order 
tomdetermine thelmicroscopic’ dissociation constants, it is 
necesSary to Obtain Gquantitative information related to 
the concentrations of the monoprotonated tautomers. A pH 
titration, the most common method of determining acid 
dissociation constants, provides information regarding 

the mean number of protons lost and thus can be used to 
determine only the sum of the concentrations of the two 
monoprotonated forms. 

A knowledge of the microscopic constants is impor- 
tant in several areas. The binding of metal ions by poly- 
dentate ligands, such as amino acids, often involves only 
one of the two monoprotonated tautomers. For example, 
the complexation of certain metal ions may involve only 
Form II (see Scheme II) or only Form III of cysteine. 
Thesbinidang by specific functional, groups is particularly 
important in peptides since the large size of many pep- 
tides prohibits the chelation of a single metal ion by 
more than a few of the multiple binding sites available. 


The ionic character, and thus chemical behavior, of 
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compounds containing more than one type of functional group 
is often affected by the degree of deprotonation of the 
individual groups. As an example, although the monopro- 
tonated forms of cysteine both have a net one negative 
charge, Form II contains three point charges whereas Form 
teTecontains only ione point charge, _ft"s likely that 

the presence Of these point charges will strongly affect 
the behavior of such a molecule, particularly at the 
10n2c staengihssetypical of biological, fluids.) Similarly, 
knowledge of the microscopic ionization scheme is 
necessary to completely describe the chemical and bio- 
logical nature of polyprotic compounds used in pharma- 


ceutical applications (14). 


A. Determination of Microscopic Acid Dissociation Constants 


Through the use of model compounds, potentiometry 
with the glass electrode has been extended to the deter- 
Mination of microscopic dissociation constants. This is 
accomplished by replacing the labile proton(s) of one of 
the simultaneously deprotonating functional groups with 
a non—-labile alkyl group, usually a methy! or ethyl 
group. The acid dissociation constant of the remaining 
functional group is then determined by a pH titration of 
the model compound and the resultant k_ is combined with 
the macroscopic constants of the parent compound to obtain 


the remaining microconstants from equations 14 through 16. 
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This approach assumes that the acidity of the group 
Litracved isenot aitected by replacing the proton (s) of 
the other functional group with an alkyl group(s). 
This assumption is the principle of Wegscheider (9) 
and its validity is dependent upon the functional group 
TOmwhich the alkyl moiety 1s attached (10,11712).. -Ccon-— 
sequently, caution must be exercised in the choice of 
model compounds. 

The microscopic constants of cysteine have been 
estimated from pH titration data by combining macroscopic 
constants with those of model compounds. Ryklan and 


Schmidt (7) assumed the pK, of S-ethylcysteine to be 


2 


equivalent to pk, of cysteine (see Scheme II) and thus 


eS 
obtained the values 8.60 and 8.67 for Pk, 3 and Pky 5 
respectively. The remaining constants obtained by these 
authors are reported in Table 1. Grafius and Nielands 
(13) used a similar approach to solve the same problem. 


These authors assumed the pK. of cysteine betaine 


2 


+ 
HSCH.CHN(CH,) 4 


COoH 


to be equivalent to pk of cysteine and the pK, of 


2 
S-methylcysteine to represent Pk, 3: The constants ob- 
tained are listed in Table 1. The model compound approach 


has been used to approximate the microscopic dissociation 


schemes of other compounds of biological interest including 
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Gincteathaone (lS \mandecitricsacid® (16.17). 

The microscopic acid dissociation constants of 
cysteine have been estimated using calorimetric data 
obtained from model compounds. Wrathall etval) (16)s noted 
tne, four kilocalorie per mole difference in the heats of 
acid dissociation of the protonated sulfhydryl and amino 
groups of mercaptoacetic acid and S-methylcysteine, 
respectively, and assumed that the AH values represented 
thesheats Ofvacid dissociation for the analogous func— 
tional groups in cysteine. This approach also assumes 
that the heat of acid dissociation, form a,qiven. functional 
group, is the same for the diprotonated form as for the 
monoprotonated tautomer and has been criticized on this 
basise (19720) se athe calorimetricyresul tsitorathe: mrero= 
constants of cysteine are presented in Table 1. A similar 
approach has been used to calculate the microscopic con- 
stants for N-methyl and N,N-dimethylethylenediamine (21) 
From caloniametrucadata:. 

More favourable techniques are those which follow 
the deprotonation of at least one of the functional groups 
dinectly, “Beneschiand Benesch (22) used the ultraviolet 
absorption of the ionized sulfhydryl group to determine 
the microscopic scheme of cysteine. Although the position 
of the absorption maximum shifted from 236-238 nm to 
230-232 nm as the pH was increased, it was assumed that 


the molar absorptivities of Forms II and IV (Scheme IT) 
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are equal s7Gonsequently, the ratio ior the absorbance at 
a given pH to the maximum absorbance at high pH represents 
ony the inaction vor "all sulfhydry leqroups which are de-— 
protonated. 
R-S_ 
ee ks pee Bae Ue = —— Ce) 
[R=SH]> + [R=-s 4] max 
Wettten gin terms of the microscopic forms in Scheme £1, 


Equation 17 becomes: 


= my, [i eer 
Sn a) Oia ek Timah Som (CNA (18) 


Equation 18, rewritten in terms of the microscopic acid 
dissociation constants, was solved by Benesch and Benesch 
using data at three different pH values to obtain Kio 
Ky3) and Ky 39° 


Several other research groups have repeated this work 


Theis results are presented in Table’ 1. 


to elucidate the microscopic acid dissociation constants 
Ormcystenwnem(S 20723) = )iCoates,) Marsdem and Riggs =(20) 
usedethe U.V. approach of Benesch and Benesch (22) ito 
StudvyetherertectsioL temperature and Wonic strength on 
the microscopic acid-base chemistry of cysteine. Their 
resultce at 25 °crand 0.1 M vonicestrength,, are resented 
pbigy diel Blre™ ORF, 

Althougmethe UV .. approachwis limited to melecules 


containing either a phenoxy or sulfhydryl group and a 
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non-chromophore as the second simultaneously deprotonating 
functional group, the microscopic acid dissociation con- 
stants of a number of compounds have been determined 
using this technique. Edsall, Martin and Hollingworth 
(11) used the ultraviolet absorption of the deprotonated 
phenoxy group to determine the microscheme of tyrosine. 
The same paper introduced the pM plot method of calcu- 
lating microconstants from U.V. absorption data. This 
graphical extrapolation technique has since been used by 
a number of authors to study a variety of compounds in- 
cluding L-3,4-dihydroxyphenylalanine (DOPA) (12), mor- 
phine.HCl (24) and phenolalkyl amines (25). Niebergall 
etpal (4) valso studiedityrosine using U.V- absorption-— 
pH data but developed both linear and nonlinear regression 
techniques to calculate the microscopic and macroscopic 
acid dissociation constants. Comparisons between the 
constants obtained by the regression techniques and those 
obtained from pM plots, for tyrosine and morphine.HCl, 
indicated that the nonlinear regression gave better pre- 
cision and lower residuals than did the other methods. 
Discussions of various calculation procedures appearing 
in the literature will be presented in Chapters III and 
IV. 

In an approach similar to that used in the U.V. 
method, Elson and Edsall (26) used Raman spectroscopy to 


determine the microconstants of cysteine. By following 
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the intensity of the S-H stretching band near 2580 cm. 
the fraction of the sulfhydryl groups which were depro- 
tonated at a given pH could be determined. Although the 
experiments were performed at relatively high ionic 
strengths (0.5-1.0 M), the results obtained compared 
favourably with those of Benesch and Benesch (see Table 

1). 

Of the preceding methods, those in which the depro- 
tonation of a specific functional group is monitored depend 
upon a particular property of the given group and there- 
fone can, bewapplied only to those molecules whach, contain 
the specified groups. For example, the simultaneous 
deprotonation. of tthe two carboxylic acid groups an qluta— 
thione or of the two amino groups in lysine cannot be 
resolved by ultraviolet spectroscopy. Consequently, a 
technique which is not dependent upon the characteristics 
of the functional groups undergoing deprotonation would 


have a wider range of applicability. 


B. The Determination of Acid Dissociation Constants By 


Nuclear Magnetic Resonance Spectroscopy 


Grunwald, Lowenstein and Meiboom (27) have shown 
that aqueous acid-base equilibria are generally rapid on 
the nuclear magnetic resonance (nmr) time scale and that 
the chemical shifts of the exchange averaged resonances 


aremluncagly maelaced to the degree of (protonation of the 
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acidic functional groups. Since all polyprotic organic 
acids have nmr active nuclei whose chemical shifts are 
dependent upon the degree of protonation, the nuclear 
magnetic resonance technique has a wide potential range 
Of Fapp lication. 

Lowenstein and Roberts (28) studied the deproton- 
ation Of citric acid using proton nmr. Since the observed 
resonances reflected the degree of protonation of more 
than one functional group, ~ a series of methyl esters was 
used to provide chemical shift data for the various 
partially deprotonated forms. However, their results 
have; been criticized by Martin» (15) "Raigler et al (29) 
studied the acid-base chemistry of the antibiotic tetra- 
cycline by nmr and Kesserling and Benet (30) used pmr 
data to elucidate the protonation scheme of a related 
compound, isochlortetracycline. Both groups of workers 
employed quaternary methyl ammonium derivatives to separate 
the chemical shift effects of the deprotonating groups. 
More recently, Walters and Leyden (31) used pmr in an 
attempt to determine the ratio of protonated sulfhydryl 
groups to protonated amino groups in monoprotonated 


cysteine. However, their results showed poor agreement 


Resonances which reflect the degree of protonation of 
amcinule une lonal «group wil be seterred LO las Unigue 
Besonances ang scnose which reiLilect the degree sor pro- 
tonation of more than one functional group will, be re-— 
ferred to as common resonances. 
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with those obtained by other spectroscopic methods and 
willbe edascussedgin detail an Chapter. 111, 

ing thesNMRy work cited thus far, the compounds 
studied have provided only common resonance data. How- 
ever, if the spectrum contains at least one unique 
resonance, the microscopic’ acid dissociation constants 
can be calculated without the use of model compounds. 

Creyi sand: co-workers (32) ‘found that the chemical 
shifts of the methyl resonances of a series of N-methyl- 
diamines reflected the deprotonation of only the substi- 
tuted amino group. By combining the fractional protona- 
tion data obtained from the chemical shifts of the methyl 
resonances with macroscopic acid dissociation constants 
obtained by potentiometry, the microscopic constants for 
both amino groups were determined. Hine and co-workers 
(33734) have dlso used. pmr chemical shift data to study 
the microscopic deprotonation schemes of a series of 
N-methyl diamines. 


The acid-base chemistry of the tripeptide 


O 
cmmmi eee a eee aes 
- 
NH, O CH. SH 
glutathione was studied by proton nmr (35). Microscopic 


acid dissociation constants for the simultaneous depro- 


tonation of the two carboxylic acid groups and for the 
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Simultaneous deprotonation of the sulfhydryl and amino 
groups were all determined from chemical shift vs pH data. 
Due to the spatial separation of the deprotonating func- 
tional groups and the rapid attenuation of proton chemical 
Sigg ieee rects st the espectrum of glutathione, contains at 
least one unique resonance for each of the simultaneously 
Geproconating functional groups. AS a result, the micro— 
scopic constants could be calculated using only the 
ehemical shitt-pH data. The macroscopic, constants, cal— 
culatedyrrom, thesmicroscopic constants via Equations 14 
and 15, agreed favourably with those in the Iiterature. 
Although the previously cited papers have utilized 
only pLOLon Nmr,—sthe, observed chemical shrits in .carbon—13 
nuclear magnetic resonance (cmr) are also dependent upon 
the degree of protonation of the acidic functional groups 
im the molecule. In proton decoupled > cmr spectra, a 
singlet is observed for each non-equivalent carbon. 
Since each carbon can be monitored, the amount of infor- 
mation is usually greater than that contained in the pmr 
spectrum of the same molecule. In addition, the range 
Of Chemical Shitts: in cmr as generally-at least=an-order 
of magnitude larger than those observed in 1 4y-nmr. Thus; 
carbon-13 magnetic resonance should also be a useful tool 
for studying the acid-base chemistry of organic molecules. 
A number of papers have appeared in the literature 


detailing the pH dependence of the carbon-13 chemical 
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Sshitt. (260-45) but microscopic acin dissociation 
constants were not determined from the data. Freedman 
et al (44) studied a number of amino acids, peptides and 
proteins but reported only the macroscopic constants of 
histidine. Fairhurst (45) attempted to determine the 
microscopic acid dissociation constants of cysteine and 
several related compounds from carbon-13 chemical shift 
data using nonlinear regression techniques. Since each 
of the carbon resonances observed reflects the deproton- 
atton Of both £unctional groups in the microscheme, 
estimates of the intermediate chemical shifts are re- 
quired. However, model compound data proved unsatisfactory 
as a source of intermediate shift values and only macro- 


scopic constants could be determined. 


Cc. Overview 


Of the techniques for the determination of micro- 
scopic acid dissociation constants discussed in this 
introduction, the model compound “approach lacks accuracy 
due sto vies inherent assumptions and the use of ultra— 
violet spectroscopy is limited to molecules in which one 
of the two simultaneously deprotonating functional groups 
is a chromophore. In addition, the calculation procedures 
presented anethe literature = endeto rely on uncertain 
estimates of vital parameters or are restricted to systems 


in which the ratios of the microscopic constants fall into 


a narrow range. 
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Since the nuclear magnetic resonance technique can 
be applied to molecules containing a variety of functional 
groups, a study was undertaken to determine the potential 
OPED re sPCCeLOsecopynitor the elucidation of acid-base 
equalgbra alton’ theimolecular level. oChaptersI11 contains 
the theoretical background and the development of methods 
fon thesdetermination of microscopic and macroscopic acid 
dissociation constants from nuclear magnetic resonance 
chemical shift data. The application of the methods 
developed» an Chapter If] is ailustrated in Chapter IV by 
tne determination of the acid dissociation constants of 
ethylenediaminemonoacetic acid and lysine; two compounds 
which are not suitable for study by the U.V. and Raman 
methods discussed previously. In addition, the data ob- 
tained for these compounds is analyzed by several calcu- 
lation methods from the literature and the results are 
compared with those obtained using the methods developed 
tnaChaptcern L111. 

Chapter V presents the results of a study of the 
determination Of acid dissociation constants for a series 
of a,w-diaminocarboxylic acids and two lysine dipeptides 
from both proton and carbon-13 nuclear magnetic resonance 
data. The acid-base chemistry of these compounds is dis- 
cussed in the light of the previously unknown microscopic 
ACidsdlSSOCLlatLon Constants. 


As part of the study of the application of CMR to 
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the elucidation of the acid-base chemistry of polyprotic 
molecules, a linear parameterization scheme was developed 
to predict the carbon-13 chemical shifts of amines, car- 
boxylic acids and amino acids in aqueous solution. The 
development of the scheme and the parameters obtained by 
the analysis of forty-five model compounds are presented 


and discussed in Chapter VI. 
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CHAPTER: i 


EXPERIMENTAL 
A. Chemicals 


i) Proton Magnetic Resonance Studies 

L-lysine monohydrochloride (Eastman Organic 
Chemicals), L-ornithine monohydrochloride (Nutritional 
Biochemicals), a-glycyl-L-lysine monohydrochloride and 
L-lysylglycine monohydrochloride (Sigma Chemical Co.) 
were all used as received. The ethylenediaminemonoacetic 
acid (EDMA) was prepared by Dr. C.A. Evans using a liter- 
ature method (46). 

A stock solution of tetramethylammonium (TMA) 
nitrate was prepared by titrating a 25% aqueous solution 
of TMA-hydroxide (Eastman) with nitric acid to a neutral 
DH. 

In the pmr titrations, the pH was adjusted with 
concentrated hydrochloric acid or potassium hydroxide 
pellets. Ionic strength adjustments were made using 
reagent grade potassium chloride. 

ii) Carbon-13 Magnetic Resonance Studies 
All chemicals were used as received from com- 
mercial sources. pH adjustments were made using concen- 


Eraoted eHelvor KOH pellets. 
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1ii1) Potentiometric Studies 
Theelysine salt was used as “received! “The titra-— 
tions were carried out using a carbonate-free potassium 
hydroxide stock solution prepared by the method described 
by Albertrand™=Serjeant (47). The KOH*solutions were 
standardized using dried potassium hydrogen phthalate 
(Matheson, Coleman, and Bell, alkimetric standard) “as 


described by Vogel (48). 


Be De Measurements 


The pH measurements were made uSing a digital pH 
meter (Orion Model 701 or 801, Fisher Model 502) equipped 
with either a standard glass electrode and a porous cer- 
amic junction, saturated calomel electrode pair or a 
microcombination electrode. The pH meter was calibrated 
using Fisher standard solutions of pH 4.00, 7.00 and 


10.00. All pH measurements were made at 25 + 1°C. 


Ce NMR Measurements 


ly NMR Measurements 


i) 

Pmr spectra were obtained on Varian A60-D or 

HA-100-15 high resolution spectrometers at a probe tem= 
Peace or 25 2 ul CG.) sBOUN spectromelers were equipped 
with Varian variable temperature control units and the 


temperature of the probe was monitored with a copper- 


constantan thermocouple. 
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Chemical shifts were measured relative to the 
central resonance of TMA and at a sweep rate of 0.1 Hz 
per second. The chemical shifts are reported relative 
to the methyl protons of DSS (sodium 2,2-dimethyl-2- 
Silapentane-5-sulfonate), the conversion being made via 


Equation 25 


Once = Stes Be Aue HS) (25) 


AIT chemivcalyshirts are reported in units Of parts per 
million (ppm), with positive shifts representing resonances 
which are downfield (less shielded) from the resonance of 
DSS. The chemical shifts are considered reproducible to 
Wei 0 HZ 6 
ii) Carbon-13 NMR Measurements 
All CMR spectra were recorded on a Varian 

HA-100-15 NMR spectrometer operating at 25.1 MHz and 
equipped with a Digilab FTS/NMR-3 system. The Fourier 
transform mode was used with broadband proton decoupling. 
Each spectrum required 2K accumulations. The system was 
focked on the deuterium resonance of Cede in a concentric 
Cavuulacyin, the 2 mn Cube. 

Chemical shifts were recorded relative to the 
Carbon resonance of 1,4-dioxane, added as internal stan- 
dard, and are reported relative Lo external IMS. ~The 
resonance Of dioxane is 67.73 ppm downtizeld from TMs. 


Positive shifts correspond to less shielding than TMS. 
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The chemical shift measurements are considered accurate 


C642 *02052ppm, 


De eOlUuCiOn Preparation 


i) PMR Studies 

The appropriate amount of amino acid or peptide 
was dissolved in distilled water to give 100 ml of solu- 
tion. Enough TMA~nitrate was added, prior to dilution, 
to make the final solution 0.01 M in TMA, which served as 
the internal chemical shift standard. Due to the lower 
concentration of EDMA studied (0.02 M), the concentration 
of TMA was reduced to 0.005 M. 

The solutions of lysine and ornithine did not re- 
Guisesany initialvpH adjustment Co obtain ar solution: of 
the fully diprotonated form. The lysylglycine solution 
was adjusted to pH 3.0 and the EDMA solution to pH 1.0 
in order to obtain a sample of the diprotonated and fully 
protonated forms respectively. All four compounds were 
titrated with KOH pellets after the initial sample had 
been taken. The glycyllysine solution was adjusted to a 
Heo Gut. sineordem tOuGbtoin dusamp le On the sully de- 
protonated sLormuprion to titration with concentrated HCl. 
With the exception of EDMA, no inert salt was added to 
control the ionic strength. Potassium chloride was added 
to, the EDMAssolutron, to bring the inatwal 1onic strength 
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All »solutsons were maintained at 25 + 21°C with a 
water bath during the pH adjustments. Samples of about 
0.5 ml were withdrawn at approximately 0.3 pH unit inter- 
Vale duringwtne course Jof the titration. 

1) CMR, Studies 

ine sEitrativons Of lysine, ornithine, 2 4- 

diaminobutyric acid and 2,3-diaminopropionic acid were 
performed in essentially the same manner as described for 
the pmr solutions with the following modifications 
(a) 1,4-dioxane was used as the chemical shift standard 
and: (b), the a2nitial, pH was reduced ,~prior to. titration 
with base, to obtain a sample of the fully protonated 
form. All solutions were 0.185 M in o,w-diaminocarboxylic 
acid. 

For the model compound chemical shift studies, 25 
EOn>0emleoG 043s Mesolution, containing approximately. 0. Lam 
dioxane were prepared. If the initial pH indicated the 
compound under study was basic, the pH was adjusted to 
aboverloewrth hOHepelletsyto Obtains the tullyedeproton— 
ated form. A sample was. taken, and the pH was then re- 
duced to a level, determined from literature pK, values, 
aires thesmmoLecuke would be fully protonated. A second 
sample was then withdrawn from the solution. For acidic 
compounds, the fully protonated sample was obtained prior 
to increasing the pH to obtain the fully deprotonated 


sample. When the compound contained more than one type 
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of functional group, the pH at which the intermediate form 
would predominate was calculated from literature pK, 
values, when available, and a sample of the appropriate 
form obtained (an intermediate form has one type of func- 
tional group fully protonated and the other fully depro- 
tonated). 

All solutions were prepared at 25 + 1°C. Approx- 
imately 2 ml of solution is required for each carbon-13 
spectrum when 12 mm tubes are used. 

iii) Potentiometric Titration Studies 

The solutions used in the potentiometric titra- 
tion studies of lysine were prepared with freshly boiled, 
double distilled water. All potentiometric determinations 
were Carried out using 50 ml (pipette) of solution and 
WecLeupebrOrmedmimider a Stream Of Nitrogen to reduce -ex-— 
posure to atmospheric carbon dioxide. For the series of 
titrations performed to determine the concentration 
dependence of the pK,'S of lysine, the pattern CLechange 
ingthe Monte Strength Of the pmr titration was duplicated 
byeadding KCl] both before and during the titratzon.) furs 
was necessary due to the wide variation in amino acid 
concentration. The titrations were carried out with 
carbonate-free potassium hydroxide at a concentration 
level at least ten times that of the lysine solution being 
CitCoaLeds 


ALL titrations were performed at 25 4 1°C. 
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Ee CaLculatvons 


The overlapping macroscopic constants of lysine 
were calculated from potentiometric data using the correc- 
tion term method of Irving and Rossotti (51); When p 
represents the average number of protonated groups per 
molecule, a plot of p vs pH for a diprotic acid will be 
symmetric about its midpoint. As a consequence, it can 


be shown that: 


PK, = PHi iq im Ory, (1.9) 


Bees e PHB ety (20) 


whereNd ~S sqreatem than zero but tess than 1 andy) is 


defined by 


eo (1+d) [Ht], _ 
y = 109 (454) 4 log fp - es | (21) 
(l-a) fH"), 4 


For the purpose of the calculation, the simultaneous 
deprotonations of the two ammonium groups of lysine may 
be treated as a diprotic acid. 

The linear least squares programme for the deter- 
mination Of macroscopic and microscopic acid dissociation 
constants from fractional deprotonation data was written 
insWATEILV -andcun "on the; IBM 360/67 ssystem at the University 


Opilberta computing: Centre. 
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The nonlinear curvefitting calculations were per- 
formed using an adapted form of 'KINET', a nonlinear 
regression programme developed by J.L. Dye and V.A. Nicely 
(52). The regression procedure used was the matrix tech- 
nique described by Wentworth (53). The programme contains 
a weighting subroutine which utilizes estimates of the 
variances of the experimental data to calculate statistical 
weights for the variables. A plotting subroutine was used 
to compare the observed and calculated data. 

In addition, KINET calculates a linear estimate of 
the standard deviation of each parameter and a multiple 
correlation coefficient. Due to the use of inexact values 
as constants in some of the calculations, the linear esti- 
mates of the standard deviations reflect the scatter in 
the data rather-than the true uncertainty of the parameter 
estimate. The multiple correlation coefficient may vary 
from zero to one, values near unity indicating a high 
degree of correlation between the parameters. 

The multiple linear regression performed to deter- 
mine the carbon-13 chemical shift parameters utilized the 
programme *MULTREGR developed by the University of Alberta 
Coppuing cer yices Mrublic: Fale Descriptions, Vol. li). 
The programme uses the abbreviated Doolittle method to com- 
pute the regression coefficients and utilizes a number of 


subroutines from *SSPLID (Scientific Subroutine Library). 
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F. Clarification of the Definition of Acid Dissociation 


Constant 


Acid dissociation constants may be expressed as 
thermodynamic constants, mixed activity-concentration con- 
Stants Or Concentration constants. To illustrate these 
three constants, consider the deprotonation of the acid, 


HA: 


HEU oe Re Da (22) 


The thermodynamic acid dissociation constant describing 


Equation 22 is defined as 


ee (a,4) (a,-) 7 ct} [AT] (Yat) (y,-) ee 
a (ain) [HA] : (Yun) 


where cor is the thermodynamic constant and is defined in 
termssOL. the activities, avs of the species in solution. 
The activity is related to concentration by the activity 
coefficient, Y;- 

When the two forms of the acid are expressed as 
concentrations, a mixed activity-concentration constant 
results. 


(aj+) [A] 


a {HAJ ae 


When the experimental determination is carried 


out under the conditions of dilute acid concentration 
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and Constant’ tonic strength, it is possible to obtain 
values of Yut from the Davies (or some similar) equation 
(49,50). The hydrogen ion concentration can then be cal-— 
culated, and the dissociation constant is reported as a 
concentration constant. 

c (H") [A] 


er airs rae 


tHe cConcertraviens of the acid and its conjugate, base in 
Equations 24 and 25 are determined from the initial acid 
concentration. 

Although the K. values obtained using Equation 
24 are dependent upon the experimental conditions, the 
calculation of mixed constants avoids the assumptions 
inherent in the activity-concentration interconversion 
and can be accurately performed under a wide range of 
experimental conditions. The acid dissociation constants 
mMeported an this thesis will be’ mixed’ activity—concentra— 
Gol COnetaltesS eam AmraCeOr Which willetactli1tate the secon 
version to concentration Constants will be given with) the 
tabulated values of the acid dissociation constants. 

In order to, calculate the concentration constants; 
it 1S Necessary to convert the pH meter readings to 
hydrogen ion concentrations. The meter readings are 


related to ayte the hydrogen ion activity. 


ph =) = 10g, 9 (444+) (2.6) 
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Thewactivity coefficient for the hydrogen ion was calcu— 


lated from the Davies equation (49,50), 


1 
“2 


log y* _ 0.512 (u) 


= 0.2) (27) 
i Webs () 


Where VYalLepresents the activity ‘coetracient Of an qion of 
charge Zs and. jl is the ionic strength, The ionic strength 


can be calculated from 
Mee ee eecean (28) 


where C; represents the concentration of ionic species i 
With charge Zs. The concentrations of the various acid 
species present were calculated from the mixed activity- 
concentration constants, determined as described in 

Chapter 111,, “The concentrationssor the inert ‘salts pre— 
sent were calculated from the amounts of acid, base, and 


salt added. 


(oS) 


_ a Od 7 r 


Aj. pete ME 30 PIE LIT TET SS ee ssnsentgat 


Inieuve. odae Ait os done ise “ace has ai i! spy ; 


nett baste fGne 


(Wag Sitet 20. "> xis doe cid. sone ne 28 
jam poo tse ate 7 =F iat Se yaneowns exif 'l je 
*OSINTIOR SPINE dapm he bakaieene oe ae 

Pci ingah in terres 7aB j syn ssicse 


Si AAP OT SAL. fi fe 43 “prj rd bap ‘bed am mi 


a zy 


abe Sen © Stha (te gonyowe wie ses bess aaa > 


CHAPTER III 


DEVELOPMENT OF METHODS FOR THE DETERMINATION OF MACROSCOPIC 
AND MICROSCOPIC] ACID DISSOCIATION CONSTANTS BY NUCLEAR 


MAGNETIC RESONANCE SPECTROSCOPY 
AS Inerodcuction 


the, Stepwise deprotonation ef asdiprotic acid ican 


be described by the macroscopic sequence: 


LS) K, 
HA => HA => aA (29) 
+ Sey ‘Naas 


Charges have been omitted for simplicity. At the molec- 
ular level, deprotonation proceeds by the two pathways 


shown in Scheme III. 


AH. 
j 
ie (LE) 
Vl 
(i) H AH, 
J FAN 
aL 


XN 


(ie) 
Scheme III 


The macroscopic acid dissociation COnstants are Trepre= 
sented by upper case K's and the subscripts 1 and 2 denote 


the macroscopic deprotonation step to which the constant 
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refers.) The lower case k’s represent the microscopic acid 
dissociation constants. The acid dissociation step to 
which a given microconstant refers is given by the sub- 
script, the last letter denoting the group undergoing de- 
protonation and any preceding letters indicating groups 
which have been previously titrated. 

The macroscopic constants are defined in terms of 


the total concentration of monoprotonated acid present. 


apt [HA] 
Ky = Deals (30) 
aaah [A] 


The microscopic constants, however, are defined on the 
molecular level in terms of the two monoprotonated tau- 
tomers. 


He 
a, +lA 4] 


Cn Sab ears ines Se) 
a) 


(33) 


= 34 
Kis “TAHT Ce 
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a,,+[A] 


aay ~ [HA] (35) 


Since [HA] represents the total concentration of both mono- 


protonated tautomers, we may write: 


[HA] = a =f [HA] (36) 


woien, when substituted into Equations 30 and Sl) leads to 
the following relationships between macroscopic and micro- 


scopic constants 


Ae = = k. + ee (Se72) 
[H. at 
i see at J 
is ? . a,+[Al a k. ak. (38) 
2 ([H,A] 1 eee og + Cee 
KK, = oe = He (39) 


For molecules in which the two functional groups are 
Ob widely dithering» acidities, only One of the two, micro— 
scopic pathways will be of importance and the macrocon- 
stants are essentially equal to the microconstants which 
describe the major pathway. However, when simultaneous 
deprotonation of the two functional groups occurs, the 
macroscopic constants describe only the mean number of 


protons lost by the molecule. 
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From the preceding discussion, it can be seen that 
in order to determine microscopic acid dissociation con- 
Stants, the acid-base equilibria must be followed on the 
molecular devel. The remainder of this chapter will des=— 
Grine the development of methods for the determinationsor 
MiCroscopic and macroscopic acid dissociation constants 
from nuclear magnetic resonance chemical shift data. The 
basis of the application of nmr to aqueous acid-base 
chemistry will be described, and the development of the 
calculation techniques will be presented. Several liter- 
aburesmethods, for the extraction of acid dissociation 


constants: from, chemical shitt data will elso be dascusseca. 


B. The Nuclear Magnetic Resonance Technique for Studying 


Acid-Base Equilibria 


In 1957, Grunwald, Lowenstein and Meiboom reported 
(27) Sthac, tn aqueous solutions of methylamine and methyl— 
amine —ydrochioride, the chemical shift of the methy. 
protons 1s Pinearly related to the ratio, of the concen= 
tration oL methylamine hydrochloride to ehe total ~concen— 
Pration Of samines sas a result, sthey concluded thatthe 
‘exchange of acidic protons is rapid on the nmr time scale 
and that the observed chemical shift values are weighted 
averages of the limiting chemical shifts of the molecular 
forms present in solution. (A limiting shift will be de- 


fined as the observed value when only one molecular form 
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is present in solution). The observed chemical shift is 


described by 


OObE = 2, 0583 I) 


where Se represents the observed chemical shift and Os 


S 
the fraction of all molecules in the form i with a limiting 
chemical shift 6. Depending upon the position of the 

nick active nucteus xelative to the functional groups in 

the molecule, a given resonance may reflect the deproton- 
ation of one or both simultaneously deprotonating groups. 
If the resonance in question reflects the deprotonation 
SeLrebethsgroups, 72 will be comnag a common resonance in 
this thesis, while those resonances which are affected 

Dy eenemOegreenOorceprotondtion Of Only one tunclional 

group will be referred to as unique resonances. 

Fractional deprotonation data can be calculated 
airectly from unique resonance data. Fractional deproton— 
aceon data can also be calculated Erom those common? reson 
ances for which the total change in chemical shift due 
torcach deprotonating functional group Can) be “determined: 
Microscopic and macroscopic acid dissociation constants 
Gon be calculated from fractional deprotonation data, bute 
only macroscopic dissociation constants can be calculated 
from common resonance data for which the individual con- 
tributions of the two deprotonating groups to the total 


chemical shift change cannot be determined. As a 
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eonsequence, <difrerent methods) of calculation have been 
developed for the determination of acid dissociation con- 
stants from fractional deprotonation data and common 


resonance data. 


Cee ther Calculacion of Acid Dissociation Constants L1rom 


Common Resonance Data 


Equation 40 can be written in terms of the molecu- 


lar forms in the microscopic scheme for the acid ee 


sobs > Og * Aro ar Ft Mggr Sarr + “re ty eee 


The subscripts refer to the molecular species as denoted 
in Scheme III. The fractions are defined in terms of the 


molecular species by the following expressions: 


a, = eae (42) 
Orr = aor (43) 
Crtr = [HA] /Cy (44) 
Ary = [A] /C, (45) 


where Co represents the analytical concentration of acid. 


Cr = [Hj AH. ] + [AH] + [HA] + [A] (46) 
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The fractions can also be written in terms of the micro- 
scopic acid dissociation constants by replacing the molec- 


ulliax forms with the expressions in Equations 32-35 


au 
i He (47) 
an+ + (kK; + ers + USES e 
hed Wee 
eae be : aeat a (48) 
+ ke. 
ant ( ‘ ze eee: a onus 
Clee Kop 
a ees 
ecient 5/5 ies abana: ai (49) 
Ht i a pan mks eee 
Kee 
Ory = mp (50) 
aut Ss (kK. + esata + ogc 


Substitution of Equations 47-50 into Equation 41 yields 
an expression for the chemical shift of a common resonance 
in terms of the microscopic acid dissociation constants 


and the Timiting Chemical Shifts of Che various molecular 


FOLMS 
a On FL aks 6 eee Car a kin Ka 0 
5 tees tee et Hees hl Hag eae, (51) 
obs ayt + (k + eye + se 


i tne slimvting eshitts of the tour molecular forms can be 
determined, the microscopic acid dissociation constant 


can be evaluated from common resonance data by nonlinear 
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Teastrsquares cunveritting to Equation 51. Since Sy and 
ory ace the limiting chemical shifttswot the fully proton— 
ated and fully deprotonated forms, respectively, S5 and 
Sry values can generally be obtained from spectra at the 
low pH and high’ pH portions of the titration curve.  How- 
ever, Ort and Csr cannot be obtained directly from the 
experimental data. When the two acidic groups deprotonate 
simultaneously, there are no pH values at which a single 
monoprotonated species is the only molecular form present. 


When the chemical shift data is analyzed by nonlinear 


kK 6 


eCurvetlLocing co Equation, 517 the: values, £0r k5Sr70 4 °xTT! 


— and (kK + eed will be strongly correlated. That is, 
errors in one parameter will be compensated for by errors 
in the other parameters. The net result will be a family 
OF SoOlmtaons;pail of which will satisty the meqression 
CracLeLia mbue WlclenG andicatlon which, 2b any, Ole the 
solutions represents the actual values of the constants. 
By replacing the microscopic constants with the 
appropriate macroscopic constants, we obtain 
+ a4Ky Say, + Ky Koo 


Z 
aro 
5 = fea eerie eer (52) 


where Ont, is @ chemical shift parameter related to the 
intermediate chemical shifts and the microscopic constants 


according to Equation 53 
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k k. 
2) foie oe ee ee Deates en ae 
Sir, (ke. tk.) °rr * free “rit (53) 
i: j sf j 


When related to the macroscopic acid dissociation scheme 
(Equation 29), Our, represents the chemical shift of the 
monoprotonated form, HA. Nonlinear regression analysis 

of the observed chemical shift data will yield estimates 
of Ki, Ky and Ont: The nonlinear programme requires 
initial estimates of the parameters being determined. 
Approximate values can be obtained from a plot of 

eee - S/S (Siy = 51) vs pH. pK, and pK, can be approx- 
imated as equal to the pH at those pH values where this 
chenicadwschift ‘ratio equals 0025 and 0.75, respectively, 


and Ont, ais 0 when the ratio equals 0.5. If the chemical 


obs 
Shitt vs pH curve contains an inflection point, and is 
asymmetric about this point, better approximations can be 


obtained as follows: 


(a) Orr, = ee at, the point of inktlect2on 
(b) pK, = pH when 9251 = 0.5 
HL I 
6 - 6 
(c) pK, ~ pH when ons __a = 0.5 
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Although the estimates so obtained may be quite crude, 
they will generally be sufficient for rapid convergence 


of the nonlinear regression programme. 


D. The Calculation of Acid Dissociation Constants from 


Fractional Deprotonation Data 


i) sSeolitery Deéprotonations 
For monoprotic acids and for polyprotic acids 

with no overlapping dissociations, any resonance observed 
is, by definition, a unique resonance. Unique resonances 
which reflect a solitary deprotonation are also observed 
WHeImeOnCUuDUNCELoOnal group Of fa pOlyprotic, acid Girters 
GLeouly einecoclirLy from the Other tunceronal Gqroupse in 
the molecule. An example is the methylene resonance of 
eheraceelc acid Ggroup,Or EDMA when) ae retlects tie sdepro— 
LOnaviwonsOr tne carboxylic acideqroup. 

In these instances, the observed chemical shift 
is the weighted average of the chemical shifts of the 


protonated and deprotonated forms. That is, 


é = a4 6 +4,6 (54) 


meee the subscripts p and d refer to the protonated and 
deprotonated forms respectively. Since only two forms 
are present, obs + om = ly which when combined with 
Equation 54 gives the fraction deprotonated in terms of 


the ‘observed shire: 
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Tee 
Wy = =o? (55) 
d Pp 


The acid dissociation constant can be defined in terms of 
fractional deprotonation data by dividing the numerator 
and denominator of Equation 30 by the total concentration 


of acid present. 


KO = —— = —— (56) 


Bewlezenagsenearelationship in equations, 55 and 56, che 
acid dissociation constant can be expressed in terms of 


chemical shift values. 


(a) 


a ene Ree ) 


obs 


When both limiting shifts can be obtained directly from 
the experimental data, a value of BS can be calculated at 
each data point for which ore falls between the values 

Or 85 and 64: However, for very strongly acidic (or 
basic) functtonal groups, one of the limiting shitts may 
be unobtainable since the functional group being monitored 
remains partially deprotonated (or protonated) beyond the 
experimental pH range for which a+ can be accurately 


determined. Alternatively, spectral overlap, due to 


other nmr active nuclei in the molecule, may hide one of 
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the limiting shift values. In either case, Equation 57 
Can be used an the appropriate Vinear form to simultan— 
eously obtain both Ko and the limiting chemical shift 


which cannot be directly measured. 


(6, - 6 ) 
st d obs 
Peed = os + poe arte K (58) 
6 = 6, + (6. - 6 Jat ” (59) 
obs d p Obs hs K. 
For example, a plot of ore vs Og - Cae 
aut 


will yield om as the iantercept and the acid dissociation 
Constant as the slope. In practice, both of the preceding 
expressions, 58 and 59, deviate from linearity as the 
observed shift approaches the known limiting shift. Con- 
sequently, if the unknown parameters are to be obtained 
by a computer-performed least squares analysis of the 
Cnenmrealwehtturcata, tes advisaeblertoyplou the wexpies— 
Sion to determine the point at which the deviation from 
linearity becomes significant. The data beyond this 
point can then be excluded from the computer calculations. 
Bradbury and Brown (54) have used Equation 59 for the 
determination of the macroscopic constants of lysine. 
Their results will be discussed in Chapter IV. 


When a large, high speed, digital computer is 
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available, nonlinear regression techniques can be used 
to extract acid dissociation constants from unique reson- 
ance data. The chemical shift-pH data can be analyzed 


by non-linear regression using Equation 60 


ad 


KikOn = ako 
HP 
obs ant t+ Se 


(60) 


This method can be used to determine the acid dissociation 
constant; alone or an combination with either or both of 
the limiting chemical shifts. However, when the limiting 
shifts are not available, corroborating calculations are 
advisable since the fraction of the total chemical shift 
range which is observable will affect the values obtained. 
Tnesapplicataon Of Bquation 60 1s 1ilustrated an Chapter 
LEAVE 
ii) Simultaneous Deprotonations 

When the system under study is a polyprotic acid 
containing two functional groups which undergo simultaneous 
deprotonation, it is necessary to express the fractional 
deprotonation data in terms of the molecular forms of an 
inci yiduategrolp rom ccheme Pill; thie fraction sCre coup 
i which is PEOLOnated, f. oo and the fraction which is 


’ 


deprotonated, £5 qr can be defined by: 
i 


| [HAH] + [H,A] 
— = 4. 

ne iii hie lAU ale otAl © ke et oe 
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ine equivalent expressions for qroup j vare: 
[HeAH. |] + [A . | 
£. = = J = + (63 
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[HA] + [A] 


Lae 65 aie is a unique resonance which reflects only the 
T 


adeqree sor Protonation of group 1, the observed chemical 


shift can be expressed in terms of f. and f. 
f 


,de 


64 = f 6 =e ak 6 (65) 


where O. % is the chemical shitt when group 4 is, proton— 
t 


ated and é the chemical shift when group i is depro- 
Ld 


d 


tonated. Referring to Scheme III, we can write: 


On = oO. = 6 (66) 


On = Oe Se 8 ok (67) 


Thus, the fraction of all molecules with acidic group i 


protonated and deprotonated can be calculated from 


id = (H|AH,] + TH|AT + TAHGT + Tal = Ory + py £64) 
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Ee ODS doplVi (68) 
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Pag ee Ee (69) 
' ey Ae 


‘can generally be obtained from the experi- 


6 and oa Tay 


Ay 


mental data as described in Section C. 


it the vobsenved chemical shitt, changes in) o- 
1,0DS 
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observed chemical shift is described by 


Dinos | eth” PA a ae a yah oe (70) 


where A. ‘ represents the total change in chemical shift 
’ 


Of 76 Cue EOnthe complete de protonazon OL Vgroup) je 


1,0bs 


The expression for the degree of deprotonation becomes: 
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ia ey eas fe Shay 


jie Sine wM ls les. ielereSlel isigrete ea is the observed shift of the 
Locally aeprocOnated Orm. 

The concentrations of the microscopic forms “in 
Equations 62 and 64 can be replaced by the appropriate 


microscopic acid dissociation constant expressions (see 


Bovataons S2-35) frormyield £. andes Inv tenis sO Late 
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microconstants. 


et eee 


Se = ao re (72) 
, 
at + ag (k; + Sa + epi 
atk. ap kk. 
Sy Ss SS ee ee 
as awe Fh aca tka t+ ke) kek 
i Hite ins yas 


The microscopic acid dissociation constants can be deter- 
mined using Equation 72 or 73 by nonlinear regression 
analysis of fractional deprotonation data. The applica- 
tO On tiesec equations, tO fractional deprotonatbionsdaca 
will be discussed in Chapter IV. 

The correlation of the errors associated with the 
estimates of the constants determined by using Equation 
72 or 73 increases as the ratio eas increases above 5 
or decreases below 0.2. However, by recalling the expres- 
Sions relating to the microscopic and macroscopic con- 
stants (Equations 37-39), the microconstants in the denom- 
iNarOor OL Equations: 7/2, and) 7/3 can be replaced by macro— 
constanus. Form Equation 7/2 wesobtain 

aytk, ap syeelss 


f. = Lip (74) 
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tiatieevaluesmror Ky and K, are known, nonlinear regres— 


sion analysis of the f. d data according to Equation 74 
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will yield estimates of k. and sae Be can then be cal- 
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We ke: 

determined vabues tor the macroscopic constants 1s not a 
limitation since the macroscopic constants can be calcu- 
lated from either a common resonance as described prev- 
iously, or a combined Se = sl data set as described 
below. The nmx spectra of most compounds will contain 
the necessary resonances. The initial estimates of the 
microscopic constants can be obtained from plots of pM, 
vs Seed where 

pM go) ee 5) 


The application of nonlinear regression analysis using 
Equation 74 is illustrated for a number of compounds in 
Chapters. lV sand “Vv. 
When both fa. d 
chemical shift data of the molecule under study, it is 
possible to calculate the macroscopic as well asthe 
microscopic constants from the fractional deprotonation 
data. The average number of protons lost by a molecule 


OLvA Gi basic acidy al a given pl, may be nepresenvced by 


h, where 


D1 
II 


({H.A] + [AH] + 2{A])/C,, 


=n + O + 2a C76) 


and fe can be determined from the 
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The a values have been previously defined (Equations 42- 
45). By recalling the definitions of f. and ¥r. given 
Eee | apyte! 


as Equations 62 and 63 we can see that, at any given pH 
Hoe) comet a + 20 = Gy + £f. GED) 


As a result, h can be determined directly from the frac- 
tional deprotonation data of two resonances. The expres- 
sion for h can be rewritten in terms of the macroscopic 


acid dissociation constants to give 


Kea + 2K_K 
— Hts © ie 
ee (78) 
ay,+ + K)a,+ + K,K 


H a TZ 
Equation 78 can be used to analyze the h data by non- 
linear regression techniques to determine overlapping 
MAaCrOSsCOpilCc acid dissociation constants. | The initial 
estimates of the macroscopic constants required by.the 
CuLcVvVelrteing programme, .Can, be Obtained from, ay plomor h 
« pH when h = 0.5 and pK 


VSmpiowinat 1s pk = pH when 


1 2 
ho= 1.5. The results obtained with chemical Shitt vs pu 


Gatagctronea vartety ot molecules ware spresented sin Chapters 


iVeand Vie une values Or Ky and K, obtained “from ehe 
analysis of h data can be used in Equation 74 to obtain 
estimates) of microscopic acid dissociation constants: 


Equation 78 has also been applied to the analysis of 


DOLenL-LoOmMetrIcmaata (55). 
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Since most nonlinear regression programmes require a 
large digital computer for their execution, a method for 
the determination of macroscopic and microscopic constants 
by the linear least squares analysis of fractional depro- 
tonation data has been developed. 

The macroscopic constants can be defined in terms 
of the molecular fractions, as was illustrated by Equations 
37 and 38. By combining these two equations, the defin- 
BELONG H h, and the mass balance of all forms of the acid, 


it Can be shown that 


ath (2-h) 


= ap aes (7,9) 
(1-h) a,.+(1-h) 


This expression is of the form y = mx + b and can be 
treated by standard least squares and/or graphical tech- 
niques to obtain values of K) and K, from h vs at Cialea. 
Data points in the proximity of h = 1 must be avoided in 
the least squares approach (51). 


ThesVvalue.oL Ka, obtained using Equation 79)cansbe 


2 
used to determine the individual molecular fractions from 
thestract onal saeprotonation data. Combining the expres= 


sions relating the macroscopic constants and the molecular 


fractions with the definition of h (Equation 77) leads to: 
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Values of Ary SO obtained can be substituted into expres- 
S1oOnselor. f. d abetsle par d (Equation’s’ 62 and 64) to obtain 

f Jr 
Orr and Ort: pance the sum of all four molecular forms 
LSe Una tyatealls pi! c,. the value. of as is readily obtained. 
They calculation) of the microscopic constants: from the 
appropriate molecilar forms: (see, Equations 32-35) as 
straightforward. The precision of the values so obtained 


LSeDeSt an thesrange (0.25. < f. i LUT 5) i eis 
lA 


eed 
method has been used to evaluate the microscopic constants 
of several molecules and the results are presented and 


discussed in Chapters IV and V. 


Eo Discussion 


Ormthe methods just -deseribed for the determindgtion 
of macroscopic acid dissociation constants, it has been 
found that nonlinear regression analysis using either 
Equations 52 01 78 will yield precise values whichsare san 
good agreement with potentiometric results. Having ob- 
tained estimates of the macroscopic’ constants, fractional 
deprotonation data may be treated, using Equation 74, in 
avsimulanimanners £0 obtain microscopic acid “dissociation 
constants. This approach will be illustrated with a 
number of compounds in Chapters IV and V. The remainder 
of this chapter will be devoted to a discussion of some 
of the methods for the determination of acid dissociation 
constants, from chemical shift data, which have appeared 
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Cohen and coworkers (36,44) have described a curve- 
fitting method for determining macroscopic constants from 
chemical shift titration curves which is based on a com- 
posite Henderson-Hasselbach equation. These authors have 
applied this method to the determination of the macro- 
scopic constants of histidine and related compounds from 


pmr and cmr data. Their model equation is of the form: 


Sea gh WOON oer ie MS EEE 1G PO oe ae ee oan 


where oe is) themchemical shitt of the fully protonated 


in 
form and PK and A. are the acid dissociation constant 

ang the total change in the chemical shitt for the deh 
PVOLOn Jeransition.. The chemical shittvtitravion curves 
are computer fitted to obtain both PK, and A. for each 
deprotonations For wa diprotic acid, when i = 1 the 
expression (LoPH-PRi) (1 + 10PH-PKi) assumes the fraction 
Ot thewacia. in the monoprotonated orm jbo, be [HA] / (HA) + 
[HA] iteembe woverlapping dissociations occur, the fractaon 
of the acid in the monoprotonated form is [HA] /([H,A] + 
[HA] + [A]), and macroscopic constants determined by 
Equation 81 will be in error. For consecutave macroscopic 
constants which vary by more than 3 pK units, (10 Boia, 
Ct Toe tee et) will represent the fraction from which the 
ith proton has been removed to within 1 ppt. However, as 


the separation between pK and PK ay decreases, the errors 
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in the fractional concentrations of the various proton= 
ated forms will increase. Figure 1 illustrates the errors 
in fractional concentrations, as defined by Cohen etjal 
for theyrirst and second proton transitions of a diprotic 
aed, asa function Of pH for different values of the 
ratio Kj) /K- 

Recalling Equation 51, one can see that inter- 
mediate chemical shift values are required in order to 
calculate microscopic constants from common resonance 
chemvcal shift data. The technique of approximating 
intermediate. chemical Shifts with the limiting chemical 
shifts of model compounds has been utilized by a number 
Ofeauenors, (26-317,33)2 for example, Rigler eu al, (29) 
and Kesselring and Benet (30) used a chemical shift of 
tetracycline methiodide to correct for common resonance 
effects in tetracycline and isochlorotetracycline respec— 
tavelly. The correction applied was small enough that 
its effect on the eres may be considered marginal, 
Particularly when compared to the errors introduced by the 
assumptiens and a2nterpolations involved in the method of 
Calculer von. (30)F: 

Martin (16) has pointed out, when discussing the 
ionization of citric acid as determined by nmr (28), that 
small errors in the values chosen as the intermediate 
chemical shifts may lead to large errors in the calculated 


microscopic dissociation constants. This point can be 
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FRACTIONAL DEPROTONATION ERROR 
oO 


Rigure 1. 


pKy-2. 0 pK- pk, so pKi#)— pKet2—oipK,#3 


Fractional deprotonation errors inherent 
inethe modelsused fora diprotic acid in 
Equatvon 81 (367,44). 
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applied to most attempts to combine model compound data 
with common resonance chemical shift data in order to 
determine microscopic acid dissociation constants. In 
general, the success of this technique is limited to sys- 
tems in which the applied correction is quite small. 

Due to the interest in determining the macroscopic 
acid cVssociation Constants of cysteine (5—8)1318,20,22, 
22,20)e ne Gecent attempt by Walters and | Leyden (31), to 
determine the relative acidities of the sulfhydryl and 
ammonium groups in cysteine deserves comment. These 
authors obtained the chemical shift titration curves for 
the carbon bonded protons of L-cysteine, S-methylcysteine 
and L-cysteine methyl ester. The change in the chemical 
ShtLt On Cysteine, Gue tO Lhe additvon, Of acid =o the 


completely deprotonated form, was expressed by 
INS neal (82) 


where f. represents the average fraction of the ith site 
which 1s protonaced and A; the appropriate protonation 
Shieteconstant. A protonacion Shiftyconstent was deLined 
eesthe met ertect om the protonation) Of a) particuals 
functional Group on the observed chemical shifts or: the 


=€H= and) —-CH— pErOvon resonances. 


2 
Thesprotonadtion shite comstant of the sulfhydryl 


Group, & for the methine proton of cysteine was 


SH’ 


calculated by subtracting the total change in the observed 
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chemical shift of the CH resonance of S-methylcysteine 
from the net change observed for the equivalent resonance 
of cysteine. Similarly, the constant for the carboxyl 
group was calculated from the data for cysteine and the 


methyl ester. Dar and “cook were then subtracted from 


the total observed change in the chemical shifts of cys- 
teine to obtain the protonation constant of the amino 
group. Similar values were obtained for the methylene 
resonance to give a second equation of the form of 
Equation 82. A third relationship defined the normaliz- 


aLLOnVOLPenest ractional “values according to 
n = tf. (83) 


where n is the number of equivalents of acid added per 
equivalent of amino acid. The estimates of f. at various 
values of n were obtained by minimizing (A6 - i no 
fon both resonances. (31). 

Careful examination of this approach raises several 
Questions regarding its application fo cysteine. ~The 
methylene protons of cysteine are magnetically non- 
equivalent at higher pH's; the observed resonance patterns 
“undergoing a rapid transition from A,X to ABC at pH 8.0. 
Sincegchic CH, group of S-methylcysteine does not exhibit 
similar effects in its pmr spectra, the assumption that 
the chemical shift parameters of the methylene resonance 


of S-methylcysteine can be used as approximations for the 
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parameters of cysteine is tenuous. Since, for deprotonated 
L-cysteine, the trans rotamer is preferred relative to the 
two gauche configurations (31), the observed chemical 

snare or the methylene resonance will bé an unequally 
welghtedwaverage of the chemical Shifts of the three 
rotameric configurations. Consequently, the protonation 
shift constants for the ~CH. - BeSOnance Of Cysteine can 
only be calculated when the populations and chemical 

Sittts Of the three configurations can be taken into 
account. 

Duey to Ene detinite separation, of the carboxy! 
deprotonation from the pH range over which the sulfhydryl 
and ammonium deprotonations occur, the use of the cysteine 
methyl ester data is unnecessary and can be seen to 
PHCELOGUCe, Crroredinto the calculation, OF the: preronatton 
sities constants... Using the spectral parameters fon 025 i 
Cysteine and S-methylcysteine given in Table I1 of 
Benerencemiol), the protonation. siith constants were cal 
culated and are presented in Table 2. The constant for 
the carboxyl! deprotonation is ‘calculated solely from the 
eysteine data below pH 5-0... "Difficulty in obtaining the 
fintewigquchemmealeshitt Of the Tully protonated Lormyot 
eylce ine Willenotwartect they calculation Of tne stractions 
of the monoprotonated forms. 

The change in the observed chemical shift of the 
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TABLE 2 


Spectral Data and Protonation Shift Constants for 


L-Cysteine and S-Methylcysteine@ 


spectral Data 


ne ae Nees as Nea Compound 
a B07 S06 473 Cysteine 
2 2292 2190 -408 
3 + LOO - 100 eae] 
e cus 
L 2b? soz S-Methylcysteine 
2 hy 2 -462 


Protonation Shift Constants®© for CH Protons of Cys 


£ 
ne 36 (325 
hee 52 (.53)F 
NHS > : 

£ 
SconH nA? (.45) 


a) Data from Reference 31, Table I1; 

b) Number of equivalents of acid added to completely 
deprotonated form. 

co) Untes*o£ ppni: 

d) Upfaeld resonance of CH5. 

e) Downfield resonance of CHp>. 

fj Erom lable IV, Reference 3. 
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groups can be described by 


A6 = Cot oH ue “ut nut (84) 


where Eon and “wut ares the fractrvons of the sulfhydryl 
and amino groups which are protonated. When the average 


number of functional groups protonated, per molecule, is 


lina ty fon + "wat 


Equation 84 yields 


= 1. Combining this expression with 


Re eeeh 
one “nat 5 “ca 


By utilizing the data in Table II of reference 31 
and Table: 2-O0f this thesis, fxH3 can be estimated to be 
0.704 and fon to 0.296, when one equivalent of acid has 
been added. The ratio, R, of protonated amino groups to 
protonated sulfhydryl groups as thus 2.4. This value 
varies considerably from the estimate of R = 0.60 calcu- 
lated by Walters and Leyden but is in good agreement with 
fhe results’ of asnumben of Obher studies (20,227,720) 2) The 
difference in the two estimates of R arises primarily 


from the inclusion of the methylene resonance data in the 


calculation by Walters and Leyden. 
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CHAPTER IV 


APPLICATION OF METHODS FOR THE DETERMINATION OF ACID DIS- 
SOCIATION CONSTANTS BY NUCLEAR MAGNETIC RESONANCE 


o PECTROSCOPY 


A InLroductk von 


In this chapter, the methods developed in Chapter 
TAIT are applied to the determination of the microscopic 
and macroscopic acid dissociation constants of ethylene- 
diaminomonoacetic acid (EDMA) and lysine (2,6-diaminehex- 
anove facta) #irom mmr chemical shire vo pu data. —Both 
molecules contain simultaneously deprotonating ammonium 
groups and cannot be studied by the other spectroscopic 
technigues described in Chapter I. 

Literature methods for the determination of acid 
dissociation constants also will be used to analyze the 
fractional “deprotonation data of EDMA and lysine: Since 
the acidity ratio, for the two ammonium groups, is about 
2 in EDMA and 8 in lysine, analysis of the fractional de- 
protonation data of these compounds will provide an indica- 
tion of the reliable range of application of the method 
used. The constants obtained by the literature methods 
will be compared to the values derived using the methods 


developed in Chapter III. 
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B. Ethylenediaminemonoacetic Acid 


i) NMR Spectra 
The chemical shifts of the carbon bonded protons 
OLVEDMA ave shown as a function Of pH in Figure 2. ‘The 


resonance due to the methylene protons of the acetic acid 


+ + 
H3 NCH.,CH.N H5CHCO5H 


group is a singlet over the entire pH range shown whereas 
the resonance pattern due to the ethylene protons varies 
fromrawsinglet at low pH to an AA"PB" multiplet in basic 
solution. “The chemical) shift of /thescentre of. the muLeip— 
fee 1s plotted in Figure 2. The change in the chemical 
shift of the methylene resonance between pH's 1 and 4 is 
due to the deprotonation of the carboxylic acid group and 
between pH's 4 and 12 to the deprotonation of the secon-— 
dary ammonium group. The chemical shitt data for the 
methylene resonance between pH's 1 and 4 is presented in 
Table 3. It has been assumed that deprotonation of the 
pEimMenyeamine Group. does not alter, .lOva Signimicantede— 
gree, the chemical shift of the methylene resonance due 
peo) the rapid attenuation of the effects of deprotonation 
as the number of bonds between the site of deprotonation 
and the observed carbon-bonded protons increases (32,35). 
The changes in the chemical shift of the ethylene resonance 
reflect the deprotonation of both amino groups. Conseg- 


uently, the methylene and ethylene resonances provide three 
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CHEMICAL SHIFT, 


Figure 2. 


NCH,CH.N 


pH dependence of the chemical shifts 
of the carbon-bonded protons of 
ethylenediaminemonoacetic acid ina 
0.02  M aqueous, solution. 
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TABLE 3 


Proton Chemical’ Shift Titration Data for the Acetic Acid 
Protons of EDMA Below pH 4.0 


Splim Chemical Shift? 
0.98 4.029 
ees 3.998 
155 3.948 
1.80 3.897 
202 3.885 
2234 3.798 
2eod 7 
2.98 sage | 
5552 oe (uk 


a) in ppm vs DSS. 
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bine 26 aar ana eit ca part 
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types of nmr chemical shift data: the methylene resonance, 
up to pH 4, reflects the deprotonation of a solitary func- 
tional group; above pH 4 the methylene resonance is a 
unigue resonance reflecting the deprotonation of one of 
two simultaneously deprotonating groups; and the ethylene 
resonance, above pH 4, is a common resonance affected by 
the two simultaneous deprotonations of the ammonium groups. 
As has been seen in Chapter III, the fractional deproton- 
ation data obtained from the methylene resonance and the 
common resonance data obtained from the ethylene resonance 
are analyzed by different methods. 
ia) Determination, of Acid Dissociation Constants 
a) Macroscopic Acid Dissociation Constants 
K) of EDMA was determined by performing a 
nonlinear curvefatiing of the data in Table 3 using 
Equation 86, which is the appropriate form of Equation 
44 
Kyi5q + ayto 5 


(86) 


De os Bart + K 


il 


Aswtnerchemical shitt of the Lully protonated form of 

EDMA was not obtained, and since the limiting chemical 
shist valuevat-sabout pH 4.0) asenot precisely defined, bork 
limiting shifts were treated as unknowns. Using the 
curvefitting programme KINET (52), pK} was determined to 
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@stimated) to be 4-075 ppm for the fully protonated form 
ana 3.709 ppm for the diprotonated form. 


Thesmacroscopic: constants, kK. and K3, of the 


2 


simultaneously deprotonating amino groups of EDMA were 
determined from the chemical shift titration data of the 
ethylene protons. Nonlinear curvefitting using Equation 
87 yielded the macroconstants given in Table 4. The 


values for 


2 
: ato y + Atk, Our =F Kj K36 54 Aas 


obs a 2 + K, + K 
Ho oot oe > 


oF and ory were obtained from Figure 2. The predicted 


Value) for Sarr, is35.002 ppm (vs DSS). initial estimates 


for Wk K. and Saar, were obtained from the chemical shift- 


Ze 3 
pH data as described in Chapter IIL. 
b) Macroscopic Acid Dissociation, Constants 
Ther Microscopic acid dissociation constants 


of EDMA were determined by nonlinear curvefitting of frac- 


tional deprotonation data using Equation 88: 


awake + ok ek 
¢ it Heel aA a ess (88) 


2,0 2 
: ayt + Kyayt + Kok, 


The results are presented in Table 4. The f. & values 
t 
were calculated from the methylene resonance titration 


data using Equation 71. A. j was assumed to be zero. The 
[? 


values of K, and K. obtained from the ethylene resonance 
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TABLE 4 


Macroscopic and Microscopic Acid Dissociation Constants 


of EpDMaAa,b 


$$ 


Constant Value Metnod or Cal culataon 
Cc d ; e 
PK, Te ee OL Equation 60 
PK, 6.085 29.0. Equation owls 
PK, Goel ee 0. Equation 52° 
eee 6.94 + .01 Equation 747 
Pk 9 : sae quation 
PK, 3 ies Oat aU! g 
PKj 53 OF SOs Ol): Equation ae 
Pk} 35 OE AO eee OZ h 
a)eeMixed=activity concentration constants. To converc co 
CONCENLLatI One Constants msubtract 09% 
D) e002) M SEDMA, us = 0-12-0016) M; 
CO) Ol agen Og (kee 
ad) Deviations are linear estimates of the standard devia- 
ticneasucal culated by KINEDT: 
e) Nonlinear curvefitting of chemical shift data. 
Ej SNonlanearecurvetitting of fractional, deprotonation daca. 
Gimme eee 82 Seer 2 
oe ee ae ee -1 
i) The subscripts denote the functional groups undergoing 


deprotonation, as previously described. The carboxylic 
acid group, the secondary amino group and the primary 
amine -cOups are: denoted as 2,92, ana 3 respectively. 
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data were supplied as constants. The MiLCLOSCOpIC acid 


dissociation scheme for EDMA is illustrated in Figure 3. 
C. Lysine 


i) NMR Spectra 

Lysine is an a,w-diaminocarboxylic acid with 
three titratable protons; As can be seen from the acid 
dissociation scheme in Figure 4, the deprotonation of the 
carboxylic acid group occurs separately, at low pH, 
whereas the two amino groups undergo simultaneous depro- 
ECOnatiOn sin thesph range 7.5°to 12-5. 80he pmr spectrum 
obtained for 0.5 M lysine HCl at approximately neutral 
pH is shown in Figure 5. The triplet centered at about 
4.25 ppm is due to the methine proton and the triplet at 
approximately 3.75 ppm is assigned to the e-methylene 
DrOLons.  Thejresonance pattern vat 3.60 ppm 15) duesto the 
methyl protons of TMA .NO, which was added as an internal 
chemical shift standard... Thermultiplet pattern over 2.00 
to 2.75 ppm is the methylene envelope for the other 
methylene protons. 

Themchemical sshitt. titration data OL the merhine 
and e-methylene resonances is presented in Figure 6. The 
chemical shift changes of the methine resonance, above 
pH 6.0, are almost entirely due to the deprotonation of 
the a-ammonium group whereas the changes observed in the 


chemical shift of the e-methylene are almost entirely due 
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Figure 6. pH dependence of the chemical shifts of 
the protons on the alpha (lower curve) 
and epsilon (upper cunve) ‘carbons: of 
lysines inwa, 0.1 Meaguecusssolutcion. 
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to the deprotonation of the omega ammonium group. For 
both resonances, it is the chemical shift of the central 
peak Of the observed triplet pattern which is plotted ain 
Figure 6. 
22) ) Determinataon of Acid Dissociation, Constants 
a) Macroscopic Acid Dissociation Constants 
The macroscopic acid dissociataon constants 
for the stepwise deprotonation of the two amino groups of 
lysine were determined using both linear and nonlinear 
methods. Fractional deprotonation data was calculated 
from the methine and e-methylene chemical shift data using 
Equation 71. Ao 3 and Az 2 were both assigned values of 
-1.7 Hz, as determined from chemical shift vs pH data for 
n-pentylamine, 2-aminohexanoic acid and 6-aminohexanoic 
acid. In all three compounds, complete deprotonation of 
the amino group results in a total change of -1.7 HZ in 
the resonance of proton(s) five bonds removed from the 
Site of déprotonation. (The negative sign indicates an 
upfield shift.) The fractional deprotonation data, 
obtained from the methine and methylene chemical shifts, 
Was COMmbined according to Equation 77 to yield an h vs 
PH data set. Computer performed, linear least squares 
analysis of the h data according to Equation 88 and non- 


linear curvefitting of the h data to 
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Equation 89 resulted in the values of pK, and PK given 
in, table 5x 
Ko ate Re K 
= 23H 
h = a (89) 
aut + atk, + KoK3 


Initial estimates of K, and K,, 


obtained from a plot of h vs pH at the points h = 1.5 and 


required by KINET, were 


h = 2.5 Fespectively. 
14) Microscopic Acid Dissociation: Constants 
Thiveee (Sets Of estimates) for the microscopic acid 
dissociation constants of the amino groups of lysine were 
calculated from the fractional deprotonation data obtained 
from the chemical shift data of methine and e-methylene 
resonances. Nonlinear curvefitting of the methine frac- 


tional deprotonation data, f was performed by KINET 


Oya 


using Equation 87. The tO a data obtained from the 


e-CH, resonance was Similarly analyzed using Equation 90) 


k ame + k k 
eee: 135132 is 


w,d 2 
ait aE a4+Kk5 + KK, 


Thiesrequired values of K, and K, were obtained from frac- 


tional deprotonation data as previously described. he 
Mierescopic constants are presented in Table 3.5 eThe PK, 3 
and PK} 30 Values given tor the methine data. and the Pkj5 


and pk Values oiven for the e-methylene data were cal= 


eZ 


culated from K, and K., as) described win Chapcrer ite 
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Table 5 also contains the microscopic acid dissoc- 
iation constants of lysine calculated by the linear least- 
Squares analysis of fractional deprotonation data. ‘The 


value of K3 obtained from the fy and fh data using 


7a To 


Equation 79 was substituted into Equation 91 to obtain 


Qiy data: 


+ £ (91) 


IV Bhyie! aed’ a 
The microconstants were then calculated as described in 
Chapter Iii. The constants reported in Table 5 are the 
median values calculated from the data points in the 
une ead 2 ay) 


The variation of the macroscopic acid dissociation 


Fange 0.25) £ 


constants of lysine with concentration has been studied 
using potentiometry. The concentration of lysine was 
varhed from approximately 0.2.to. .005.M..) The sonic 
Strength was controlled using KC] so that the changes in 
PoOuLeCIsErengem Cunning the course Of tie tLtrallom were 
Similar to those observed in the chemical shite Crerataons 
The macroscopic constants were calculated from the pH 
data using the correction term method described by Irving 
and Rossotti (51). The results are presented in Table 
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TABLE 6 


Vardol1on il -enesMacroscopic Acid Dissociation Constants 
Orubysine With Concent rataon 


Concentyation ikon! le) sie Dic 
(Molar) Serengtuh s(M) > LS he ; 
~194 pe OS 5 SS: 105380 
UST e223) Oneal LORS 
-049 G5 = RS) Oe | LOL SL 
s0L0 a2 oe Ne Paak 10.68 
.005 see 26 Tag ail ORO 
OOS oi OTs). 10.60 


a) Calculated using Equation 28. 

DipeMixcdsaceivity Concentration, constants. 

ec) Calculated using the correction term method described 
bDyeuLVvandegandsROSsote. (515) 
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De _ Discussion 


No literature values are available for the macro- 
scopic or microscopic acid dissociation constants of EDMA 
Or for the microscopic constants of lysine. A number of 
workers have used pH titration data to determine pK, and 
pK. of lysine and the published results are presented in 
Table 7. ‘Comparison of the results for the macroconstants 
presented in Table 5 with those in Table 7 shows a con- 
siderable discrepancy between the two sets of macroscopic 
constants. However, the results presented in Table 6 
indicate that the differences are primarily due to the 
higher concentrations used in the nmr titration. 

Although the application of the nmr technique is 
intended as a supplement to potentiometry with a glass 
electrode, it iS inevitable that comparisons will be made 
between the two techniques. Examination of the equations 
developed in Chapter III and applied in this chapter will 
show that the nmr method, for the calculation of acidity 
constants, does not require that the concentration of the 
acid in solution be known. Thus, the purity of the solute 
need not be known. This is of particular advantage when 
dealing with biological compounds such as amino acids and 
peptides since the crystalline forms of these compounds 
are often difficult to purify and may contain non- 
stoichiometric waters of hydration and unknown concentra-— 


tions of salt. Similarly, the presence of extraneous 
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TABLE 7 


Macroscopic Acid Dissociation Constants of Lysine 
Determined From Potentiometric, Data 


PK, pK, Reference 

9.19” 10.62» 450 

Se vieke 10.58 a2 

pene ore on Ag 

ibe eae This thesis 

ajo Mixed activity—concentration constants; 

byes 0, 00S Meivysiner HCl qi s=" 02 Me 

G) Lysine concentration not given, Ww = 021M: 

d) Calculated from the concentration constant 
usang ¥y <=) 21/715. 

e) Calculated from the thermodynamic constant 
using the value of y = 2.775 given by the 
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Sources Of acid Or base, such as the absorption of carbon 
dioxide, have no effect on the nmr calculations. 

Aside from the complex and expensive equipment in- 
volved, and the time required to perform the data collec- 
tion and calculations, the nmr method suffers from sev- 
eral disadvantages. Chemical shift measurements cannot 
be made with the same precision as is obtained for a 
Carefully peritormed potentiometric titration. Secondly 
unless Fourier transform techniques are available, the 
concentration of compound, and thus the total ionic 
strength, are higher than those required for potentiometry 
and may be a limiting factor for compounds of low solu- 
Dmieviey. 

As previously indicated, the microscopic acid dis- 
sociation constants of EDMA and lysine have not been prev- 
iously determined and thus the results in Tables 4 and 5 
cannot be compared to literature values. However, the 
a-methine and e-methylene resonances of lysine are essen- 
tially two ditterent data sets and as a consequence provide 
aneinternal check regarding the reliability of the values 
obtained for the microscopic acid dissociation constants: 
The excellent agreement between the nonlinear calculations 
using the two fi data sets and between the linear and 
nonlinear calculations (Table 5) serves as an indication 


Of the reliability of the two methods. 


Bradbury and Brown (54) used pmr chemical shift 
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data to determine the macroscopic acid dissociation con- 
stants of lysine and several methylated derivatives of 
lysine in Do0. They assumed that the methine and e- 
methylene resonances of lysine were unique, which is 
essentially correct, but also assumed that the dissoci- 
ation of one protonated ammonium group would not affect 
the acidity of the second protonated ammonium group. 
Examination Of the wesults in Table: 5 aindicates that tthe 
second assumption 1s not valid. Consequently, Bradbury 
and Browns use of Equation 59 to calculate the  macro-— 
scopic acid dissociation constants is highly questionable. 
The pM method has been the most frequently employed 
technique for the determination of microscopic acid dis- 
sociation constants from fractional deprotonation data, 
DaLcicularly that obtained by wltrayvrollel spectroscopy. 
This method was first presented by Edsall, Martin and 


Hollingworth (11) and is based on the function M.- 


Mee eee Cel eh ae (92) 


Examination of Equation 92 indicates that, when 


ayOLOt. ot pM, vs f. 1S extrapolated «tO Zero, pM, = 
Li 


d 


pk and when extrapolated to one, pM; = oe The 


Values. Obtained £or k; and eo) can be treated as known 


quantities in order to determine the remaining constants. 
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ida plot yor pM, vs See at the point where Ee =O Oe 
pM; = pi. (AS sa, result, when f. am 0.5 Equation 92 can 
v 
be rearranged to yield: 
aa at (kK, ~ ay+) 
; oe es C33)) 


H+ 44 


In the instance where pk. or ea is approximately equal 
Loy the pH at ed = 0.5, Equation 93 becomes dependent 
upon the determination of a relatively small difference 
between two large numbers and consequently suffers a loss 
in precision. As an alternative, the value of ea fOr 
which pM, = (pk, + Shes may be substituted into Equa- 


C7Onwl2 co, yield: 
ke, bee logn = = ek (94 
Pe ae aa een esi ) 


Both of the preceding equations are essentially 
Ones point, Calculations and are Chus subject tor the same 
errors and criticisms (51) as the n = 0.5 technique used 
for determining formation constants. ~As has been prey— 
1Ously moted (11,57) the extrapolation to determine k and 
keueas accurate, tor only a lamited mange of valucs tor the 
raAcLO ae When ise >> key thewextrapolLarlonyopsene pM, 
vs aed plot sto. obtain ei is Straightforward out yehe 
steepness of the function as fiid approaches zero makes 
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Be are key pk can be determined with confidence but the 


rapid change in the slope of the curve as f. approaches 


7a 
unity makes an accurate evaluation of Pk. ALi hvewdce 
Edsallsetval (11) determined that for both ‘constants to 
be evaluated with a reasonable degree of accuracy, the 
ratio ss must lie between the values of 5 and 0.2. 
Tne values fOr the microscopic acid dissociacion 
constants of EDMA and lysine have been determined from 
pM; plots and are presented in Table 8. For lysine the 
extrapolation and calculation was performed for the frac- 
tional deprotonation data obtained from both the a-methine 
resonance and the e-methylene resonance. The pM, vs Sea 
plots for EDMA and lysine are presented in Figures 7-9. 
The ki 2/13 FatioLor EPMA» 16 29 -and comparison 
of the microscopic constants in Tables 4 and 8 shows that 
the constants obtained from the pM plot are in good agree- 
ment with those obtained by the methods used in Section B 
of this chapter. A comparison of the results for lysine, 
which has a Kj 9/13 Eatio Of Jeo, ,1ndscates pool internal 
agreement for the constants determined from the pM, vs 
Bae plpte (Tables 5 and 8)=) “A possible “explanation ior 
the discrepancies observed lies in recalling exactly 
what information is present in the nmr chemical shift 
data. The expression for the observed chemical shift of 


a unique resonance can be written in terms of the fractvon 


of the acid which is present in the various molecular 
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TABLE 8 


iNew MaccOoscopremand Microscopie Acid) Dissociation 
Constants® of EDMA and Lysine.HCl Determined From 


pM Plots 
EDMA LYSINE 
pk 5 6.91 91340 3 Be* 
see tee 10.38 10.15 
eee 9.96 10.99 10.81 
Pky 35 On55 9.95 10,22 
pK, Blah 9.30 9.46 
pK, 10.10 10.99 10.91 


alee Muxed activity concentration, constants. 
bye eealcullated using £5 54 datasot lysine. 
c) Calculated using 2 el data or lysine: 


pM, 
[e-e) 


Piguses 7. 
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PM> vs f2,q, the fractional deprotonation 
of the secondary ammonium group, for EDMA. 
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Figure 8. pM5 vs fj gq, the fractional deprotonation of the 
alpha ammonium group, for lysine. 
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Figure 9. pM3 vs £3" ds the fractional deprotonation of the 
epsilon ammonium group, for lysine. 


ry Da 
yea. 4 
ane Sale ees es 
- ers 
‘on 
NO 
sibs 
£0 


yt | 
iZ) % 
T 
wt sme 
iio 
bet 
PS 
ain 
‘ie 
ot _ 
ret 
ca 
aes 


Formss 


ea pe ee! SBP 1s (ee) Nelgs (95) 


The magnitude of the ratio ye is indicative of the 
Prect10n, Of “all molecules which will dissociate. by a 
given pathway and thus represents the ratio Ory/%p7r° 

For a Kk ratio close to unity, approximately equal 
concentrations of the two monoprotonated tautomers will be 
present at a given pH. As the ratio moves away from unity, 
one or another of the two monoprotonated tautomers will 
begin to dominate in solution. For instance, at a Sey Ae 
Betio Of 5:1, eighty—-three percent of “the acid dissociation 
will occur via the upper pathway in Scheme III. If the 


Value, OF is always much less than that of a the 


dered 


sensitivity of the observed chemical shift to changes in 


i. 


the amount of Form III present in solution will be low. 
As a result, the accuracy with which Ky and 5 can be 
determined will also be low. 

Tn addition, the extrapolation Ob the pM plot tends 
to weight the data taken at the two ends of the titration 
Curve more heavily than that obtained in the middle of 
the titration.) in those areas Of the experiment. where 
the values of f. are close £0 either Zero Or One; thie 
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increases. Consequently, the extrapolation depends on the 
Teast. precise data. This 15 also true for BA Al values ob- 
tained from other spectroscopic techniques such as ultra- 
violet abserption. In contrast, the nonlinear curve- 
PIcCC MG techinagucwutilizes alls or etne, experimental data 
and weights the data points obtained from the different 
Pats Of the titration Curve: appropriately. 

Fung and Chen (57) rearranged the equations used 
in calculating pM to obtain a linear form which would re- 
duce the extrapolation errors for unfavourable values of 
the gt ratio. When Kk. LS elartgcmecian Bug the “extra=— 
polation to obtain ke US Stralqherorward) but, the value 


obtained for 5 is imprecise. Equation 92 can be re- 


arranged to obtain the linear form: 


at (k; = M.) = seas Se Pad (96) 


Usingechnesextrapolatedsyadue OF, k, as a known quantity, 
=7Me hepa Aa 

a plot lot at (ks M.) vs M; shouldeGqivewasstralgiitanic 

On ES lope Mg and intercept - ype ak in practice, Noweviex, 


the range of experimental values for which a plot of 


Equation 96 is linear is somewhat limited and the plots 


are very sensitive to the value of ke weed) (57) .) Plocs 
he ant (ky 5 - M,) vs M. for EDMA and of at (kj 5 = M,) vs 
M, and a+ (ky 3 oS M,) vs M3 for lysine are presented in 


Figures 10-12. The values obtained for the microscopic 


constants are in Table 9. As can be seen by examining 


we ‘ at 7 mee 


lad oil LF ETIO, 9s 
‘ae 4 i . ; vt : 
oS on ft) P : : “ . 
A po a6! sib \ 
e ay a 
r y. 5 7196 Th ta > arias a OUP) 7 
re : ' YY 7 
boris wot tg 398as Fes 
wer 4 ay a 6 1 (ae 
7 , 
awe Pr + | 
Ps ¢ 4 AOD2 = phSs 4 
2 Me 
140 a sib 
} 
eo: 
19° vr | Ae SEO ls a 
Ss 
eo 7 
<n joss ve 
b ‘ if t-te : a 
¥ yy 
eed es AT ge 
5 
i ¢ a-yrt t. Liens. 
om ' ; 
7 I 
jee ain I 
a k ’ 
Goh tare 
BiGy : 
af ‘ rey ' yA ree | hey + 
‘ : | 
| (2.4 = 
ee 
H 
v, , ’ 
| 
x i cy bas 
"yl E tial fi sf ay us _ AZ 
4 : a : 


{ daitieng nl «pda + sqmemini Blew 
| . at 


ty epee 
he gaan) she Tuas (sao: ste nkse = 
. ee 


= 
4. tq “ti 2 rad £44 #5 t esi vate aio sent 

4: | oa 

“Jectot y cial. geal i ste, ay inh whet Soll . 


‘Sy Ut = > ae oy 2°28 On 
_ 4 ~ aL ; ; ; 
7 — _ 

a  Detn2eeag ath sat) 7h 
7 a 

a 
iL 4 hi a ai: - ar 
) 

nes er - a S252 


—~ 


in 


ms | 


oa 


16.0 | 
14.0 


12.0 


17 
at (k,,-M) x 10 
fo) 
oO 


0 2.0 4.0 6.0 
9 
M x 10 
Figure 10. ayt (ki 2-M2) vs M, for EDMA. The 


slope yields ki3 and the intercept 
Sai tes don oeahe 


89 


—— nO 


efhT 


wht \ 
\ 
Sa et 
ee: 
; Z ‘Ot ~ M ae ; 
Woe mee oe “a 7 io. eRe 


qaspisgna h4- 


a 


Pe 


ts) 


1) eee) 
: eo 
€ , os i 
> 
ee ee D 
a rS Glee 7 
- : perme 


Bera Gorn 


* pele x ~ ld 


20F. 
ce) 
Qe 0 
2 
x P 
- 
= 0 
8 
ay 
_ e 
+ . ? 
i 0 
8 
-1.0 7 | | 
Oo 2.0 4.0 
10 
M. * 10 


Pigure sli vayt (kis -Ma) evs Meetor lysine. the 
slope yleldsek;s) andetherantercept 


- ki3kKi23. 


90 


oe 


\ . 
i \ 
\ ¢ | 
; [ 
7 ; * ' 
iJ we 
| “oer “- 
‘. 
Pes | i 
uh 4 


ore eal sat)’ «<M av { gh : 
exfogii -ama-bug \.% .ojong ie 
a : . 


Paqure ii. 


1.0 2.0 3.0 
M,x10" 


4.0 


ant (kag Ms) VvSeM. “fore lycime seurhe 
slope yields k,;»2 and the intercept 
= Kyok)32- 


on 


= 


a 
a. 
a] 


att Sh Layl A'h Mi BY. Aa axle re ig 
~t mic Bab aya chkoky au = 


a 


TABLE 9 


Microscopic and Macroscopic Acid Dissociation Constants® 
of EDMA and Lysine Calculated Using the Linear 
Modification of pMP plots 


EDMA LYSINE 
Pky 5 6.91 934° 9.404 
eae 7.39 10e38 10.19 
Bie 9.40 10.98 10.88 
pia 9.88 9.94 10.02 
pK, 6278 9.30 Gus 
pK, 10.00 10.92 10.88 


a) Mixed activity-concentration constants. 
b) Reference 57. 

ej Calculaced from £74 data of lysine. 

d) Calculated from £35d data of lysine. 


see 


the results in Table 9, the linear modification does not 
provide a significant improvement over the values obtained 
from the pM plots for EDMA and lysine. 

Niebergall, Schnaare, and Sugita have pointed out 
(14) that both the pM plot and the linear modification 
are susceptible to investigator bias. These authors Wave 
presented both linear least squares and nonlinear regres- 
sion techniques for the determination of microscopic acid 
dissociation constants from the fractional deprotonation 
data obtained by ultraviolet absorption spectroscopy (14). 
The microscopic constants of D,L-tyrosine and morphine 
HCl were determined by performing a nonlinear regression 


analysis of f. data using Equation 97. 


“ie! 


ce = (97) 


Equation 97 can be obtained by substituting the relation=— 
Ship in Equation 39 into the enor of Equation 74. 
Niebergali@ec al (14) obtained) values for ky Ky 
and K, kK, using Equation 97 and the remaining microcon- 
stants were determined using Equations 37-39. The fod 
data for EDMA and Aya and fd date wror Jyvsine were 
analyzed by nonlinear curvefitting using Equation 97 and 
the results are an Table 102 “the results for EDMAVare 


essentially identical to those calculated using Equation 


Jau(see Table 4)e.) Snalysis of the lysine epsilon— 
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TABLE 10 


Microscopic and Macroscopic Acid) Dissociation Constants” 
of EDMA and Lysine Determined from Nonlinear Curvefitting 


EDMA LYSINE 
Pk 6.94 9.40° 
Pky 3 d38 10.26 
Pk1 93 che (ots Loe 0 
ce 9.40 9.84 
pK, 6.80 9.34 
pK, 9.98 ORG 


a) “Mixed activity—concentration constants. 
bP Using Equation, 97. 
Cc) Erom £ data. 
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methylene data yielded results in good agreement with 
those in Table 5 but the alpha-methine data gave a nega- 
tive value for K,. In addition, the statistical 1ntorna— 
tion provided by KINET indicated higher linear estimates 
of the standard deviations of the parameters and a high 
degree of correlation between the estimates of the three 
parameters. Analysis of ced data obtained from computer 
simulated chemical shift titration curves indicated that 
HOMAInecarICunbveritcing USing Hovation 97 1S mosterelwanlec 
for systems with a low value for the ratio ayia Al- 
though quantitative limits are somewhat difficult to de- 
fine, the degree of correlation approached unity as the 
es ratio became progressively greater than 5. Asa 
result, Equation 97 should only be used when other tech- 
niques indicate a low ratio for the microscopic constants. 
Shrager~ et al) (50) have described a -computerscunye— 
fitting method which utilizes nmr chemical shift titration 
Gata to obtain both macroscopic and microscopic acid dis— 
sociation constants and have used their programme to 
determine the acid dissociation constants for L-histidine 
and several related compounds. The programme is designed 
tO curvefit data trom both solitary and overlapping dis— 
sociations: calculating macroscopic constants in the 
former cases and microscopic constants, from which the 
macroscopic constants are calculated, in the latter cases. 


Their method is based upon a rearranged form of 
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Equation 51 and is written in terms of proton association 
or formation constants rather than dissociation constants. 
For systems in which common resonances are observed, the 
intermediate chemical shift values are defined in terms 

of a fraction, p, of the total observed change in chem- 
ical) shift. As’ a-result, the microscopic acid dissocia-— 
tion constants determined by this method will be dependent 
upon the values chosen for p. This is confirmed by 
examination of the constants presented for L-histidine 
methyl ester in Table 4 of reference 56. 

There are several differences between the curve- 
fitting procedure presented by Schrager ec al (56) for 
LNewAetermination Or Microscopic Constants and tne: pro— 
cedure used in this chapter. The most notable differences 
arewthat Schrager et al vcurverit ithe observed chemical 


shift data, rather than f. data, and more important, 
J 


d 
attempt to determine both the macroscopic and microscopic 
acid dissociation constants from the chemical shift 
titration data of a single resonance. Analysis of com- 
puter simulated titration data, using equation 51) indic— 
ates that if the observed resonance is strongly affected 
by the degree of protonation of both groups, the micro- 
constants obtained will be very sensitive to the values 
chosen for the intermediate chemical shifts. When the 


observed resonance is unique or is affected almost en- 


tirely by the deprotonation of only one of the two 


Cipeer | 
3 _ ) Ce 


-s a 
rs Za Poe a . ra a 
Abpea lsc-s2 = UeIeTg ig emits | 
ray , © 7 


’ 
hs me ca 
? how ssecdiyvo 28 = 
} 
4 Tu, Se Py cys 
—arTigkl’ DEheaao 
* ¢ - 
a i , m 44 , 
[ aa RAS 
f » iw 
+> 
Te : : 
— 1 
} J 
2 iii 
' ' 
poy st 
. cane 
~ + 


¢ : y Ron . ho 
STA al TA ae =i Jc ae talnesad 
J ; | ; i a = — | 
“Toe J slerta .oonakomemalyhis « to 238b Godenares 


7 
a i . - , 4 
~afonr fC, (Oar RusS Bais <j 
; 
a _— 
3 ~TD Palit) re af 
az) ee rit S 1 dish 


: aa _. a = ” 7 > 
ht iy 2qres rod aD, paar 


Fi - _ 
1 é ’ 4 - oF, > 
aay Lay ett oF, avi2be: ‘botit ¢FdO PRE 


Se 


aiid maaW - edt! oe 


: i 


>) sa si ead 9 hose 


is ca 
- ‘io A a 
"a 1 vt 2s, in _ p S 


deprotonating functional groups, the reliability of the 
Gstamates Obtained for =the microscopic constants will) de— 
Grease ac) tne ratio ba increases. Consequently). non 
linear Curvetitring of chemical shift data to Equation of 
can be expected to yield reliable estimates for the macro-— 
SCopic and microscopic constants of only a limited number 


of molecular systems. 
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CHAPTER V 


THE DETERMINATION OF MACROSCOPIC AND MICROSCOPIC ACID 
DISSOCIATION CONSTANTS OF a,w-DIAMINOCARBOXYLIC ACIDS 


AND LYSINE DIPEPTIDES 


Ko Introduction 


ine acid-base chemistry cf the majority or amano 
acids and peptides has not been characterized on the 
molecular level. in this chapter, the results of the 
determination by nmr of the macroscopic and microscopic 
acid dissociation constants of a series of a,w-diamino- 
carboxylic acids and of two lysine dipeptides are presented 
and discussed. As can be seen from the microscopic acid 
dissociation scheme in Figure 13, the compounds studied 


all contain two simultaneously deprotonating amino groups. 


Be Results 


i) a,w-Diaminocarboxylic Acids 
The chemical shift titration curves tor the 
a-methine proton resonance and the y-methylene protons 
resonance of ornithine are presented in Figure 14. The 
data for lysine has been presented previously (Figure 6). 
The chemical Shitt of the central peak of the observed 


triplet is plotted in units of ppm vs DSS. The chemical 
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Figure 14. 


pH 


PH dependence of the chemical shifts 


Of the protons onvthe alpha jandydelta 
Cacoons CL Ornithine, ©The formula or 
ornithine is H5NCH5CH,CH.CH (NH,)CO,H. 
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LOA: 


shift changes observed for each of these resonances are 
assumed to be almost entirely due to the deprotonation of 
one functional group and the microscopic and macroscopic 
acid dissociation constants have been calculated from 
fractional deprotonation data using the methods described 
forelysine an Chapter iV. “The results for lysine and 
ornithine are presented in Table 11. 

the  CMRechéemical shift data, as ppm vs antéernal 
avoxane for 2,3-diaminopropionic acid; 2,4-diamanobutyric 
acid, ornithine and lysine is shown in Figures 15-18. As 
is indicated by the data in Chapter VI, the effect of the 
deprotonation of an ammonium group will be reflected in 
chemical shift changes in the resonance of a carbon five 
bonds removed. As a consequence, the CMR chemical shift 
data of the alkyl carbons was treated by the common 
resonance techniques, described in Chapters III and IV, 
foretnevdelerminatlon of MaACEOSCOD1C constants. | [newre— 
sults are presented in Table 12. 

Attempts to calculate the microscopic acid dissoc— 
jation comstants from carbon—-l3 chemical shift data using 
Equation 51 and intermediate chemical shifts predicted 
from the linear parameterization scheme in Chapter Vi 
were unsuccessful due to the sensitivity of the calcula- 
tion to the uncertainties in the predicted chemical shifts. 
However, the microscopic constants of 2,3-diaminopropionic 


acid, ornithine and lysine have been estimated from 
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TABLE ll 


Macroscopic and Microscopic Acid Dissociation Constants 
GOfVOrnI thine. HCl and lysine sic] 


ORNITHINE> LYSINE© 
a 
pK, 8.85" 90277 
a 
pK, LOR OL LOW Se 
£ 
pk, 5 emotes) ceno1e Secu le) | eeyncol 
Pk, 3 aU eerf) LO 33. 107530 
Pk} 93 LOS O76 10/62 = 1 Ongey 
Pk} 30 he 97268 9.74 SA) 
a) Mixed activity-concentration constants. To convert 
ie, approximate concentration Constants, subtract ~09:. 
Dime UO Meornathone Cle. i802) 0s —=. 2420. 
CU oonM ysimesnCly = 0.86 -— 37/2 Me 


d) Nonlinear curvefitting of fractional deprotonation 
data to Equation 78. 

e) The linear estimate of the standard deviation was 
e002 Ofsa pK Or pk unt ror ald sconstants) listed. 

f£) Subscript 1 denotes the carboxylic acid group, 2 the 
a-amino group and 3 the w-amino group. 

g) Nonlinear Curve fitting OEaES /d data to Equation 74. 


Porson chine iy cin a err) Hz, for lysine pri 
= WS Th Asa 

h) Nonlinear curvefitting of f3,q data to Equation 90. 
For ornithine oF aie Oe 0Niz efor ely ine A3 2 = 
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CHEMICAL SHIFT, ppm vs DIOXANE 


Figure 16. 


pH 


pH dependence of the chemical 
shifts of the alkyl carbons of 
2,4-diaminobutyric acid in a 
0.185 M aqueous solution. 
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pH dependence of the chemical 
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ornithine in a 0.185 M aqueous 
solution. 


105 


-40 


1 ’ t 
w Ww ww 
& O oe) 


' 
(eS) 
NO 


CHEMICAL SHIFT, ppm vs DIOXANE 
N 


Ze) 


Figune, 16. 


6 8 10 12 
pH 


pH dependence of the chemical 
shitts “of the alkyl icarbons of 
lysine in a 0.185 M aqueous 
SOLUtTON:. 


106 


Tee | 
seul ine 


vy yi hy 


107 


AdFAzIDV “aVd-y'Z 103 WH Sy* = rt pus skT pue ‘uio ‘ayq-¢'7 103 WW ¢° 


S0°F88°97- 


90° FTE °6E- 


90° F£9° Hy- 


€0°+86°ZE€- 


$0°+08°OT- 


T8°OT 


T8°OT 


T6°OT 


90°+£0°TT 


TO" FT6°OT 


OT’ ¥66°0T 


70°+68°OT 


ST°+62°0T 


12°6 


Te"6 


07°F0S°6 


LO°F10°6 


80°+7E°6 


T0° ¥62°6 


S0°+97°6 


* (2S) 


ud yoea worz 60° 39819qGNs ‘sqUBISUOD *>dU0D 


= 09°OT 68°s 

= 89°0T €8°s8 

= S$9°OT 98°8 
£0°F98°9Z7= €0°F69'OT L0'F88'R 
80°F0L°7H- = L0°FS9°0T LO°FhZ°g 
€O°F9E°SE- 40°F8S°OT T0'F9g"°g 
LO°¥60°TTI- 8T°FS9°0IT 0°F76°8 

Vi fy yd 

ESS TE RATE) 


IPC SIFFUS TEFWaYD ET-UOqIe) WoIy paze[Na[e) SPfov 


*(0$*67) uotaenba sazaeg ay3 YITA paqetnoq{es 8 UaTIFJjaoo 
= M1 103 ‘49 *yay Jo squeysuos dFweudpowIsy ayy Woz *oTeI SqUBISUOD pax, 
*“Jue{8U0D ay UPYZ 133¥a13 SBM JULysUOD p|zzTJ |ayI JO IoIIAa piepueqs 


*9UexOFp [TeulajuyF sa uadd 


UOFIBFASP PlBpuvzs ay. Jo ajewyysa aeauT{ ayy sf pajaodai AQuyeziaoun sayz 


= Se°OT 92°83 
= 97° 0T T7°8 
fO°F86°8Z- =BI°FIS*OT 0° FZH'g 
VO°F9E*HE- TO°FZP°OL 7zO°FO"°S 
= P P 
MEE fod yd 


Poy oFak3nqouruergq—y"z 


cl alavi 


*xoidde 03 31aau09 ol 


DFTAXOqivooupuREFg-m'D jo po 2ueasuo 


voFJepoosstTq pyoy ostdoosoasey 


“SQUEJSUOD UOFZBIQUaIUOI~AQTATIOe paxTy 


= €5°6 LL°9 
= €S°6 8L°9 
ZO"+LO"HZ-  TO"FES*6 = ZO" FLL'9 
€O°0Z"ZT- ZO" FES*G = 20" FELD 
Vu € z 
ae xd a 


pfoy oyuoydoidouymerg-¢"z 


> 
_ 


stan a ai shee 
ae hs 


7 i a 


: - — = aeas 5 
- - ane - 7 
7 . : &_43: —— = ee ee 
a) a, oe , 
= 7 : —z, —_— = © 


Rt .bi oY Le : 7 


: a pas _ 7 
a id o : es 
es ee ve Sewie ‘amt, CS ae . - 


oor 1. bree = 2 ; : 
ee T° a re 4 — 14 Te. & : a of i. 2 
: &e,! 7 
7 ~ a if is 
- 
_ 
_ = 
=. - — : te. i ; os 
wae 
+ 
on = Sore vs fal 
~— wa vate ,ou7 5 c : 
7 a et : " > : oe i + ise ats -eF a 
> _ - 
1 ee 
> » oar ts 
_ : : oa 1,49 wt 
7 - = . 
ia be pris: GAGs F307 Kes, = bre Ges bin i oat ed 
2 | 7 
a 
» 


TOL 


108 


fractional deprotonation data using ge values predicted 
from the model compound data in Chapter VI. The predicted 
ne Values and the results are presented in Table 13. 
ii) Lysine Dipeptides 
The pmr spectra of a-glycyllysine and lysylgly- 
cine contain an essentially unigue resonance for each of 


the two simultaneously deprotonating amino groups. For 


e-glycyllysine the chemical shifts of the glycy! methylene 


+ oF 
H.NCH,CONHCHCO.R aries 
(CH,) a We ) 3 
i (i 

+ + 
NH3 NH} 
a-glycyllysine lysylglycine 


and of the e-methylene of the lysine residue are presented 
in higure =) 9. ee inewglycyl: protons Give rise, to a singlet 
over the entire pH range in which the ammonium groups are 
deprotonating whereas the e-methylene of the lysine frag- 
ment yields a triplet which exhibits some broadening at 
the higher pH values. Both the a-methine and the 
e-methylene of the lysine fragment of lysylglycine give 
rise to a triplet pattern over the pH range Of Interest 
and the chemical shift of the central peak of each triplet 
is presented as a function of pH in Figure 20. These 


chemical shift titration curves have been analyzed by the 
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TABLE 13 


Microscopic Acid Dissociation Constants of 2,3-Diaminopropionic 


Acid, Ornithine and Lysine from Carbon 13 Chemical Shift Data.@ 


ane ORNS Lys? 
. | ce ce Ce Ce Cy 
eae 7068 7.07 8.93 8.92 9.35 9.30 
PK, 3 76) 7.09 9.63 9.67 10.02 10.43 
Pie 9.25 9.24 10.52 10.53 eT 10.74 
Bhi 9.21 9.21 9.82 9.79 10.06 9.66 


a) Mixed acidity-concentration constants. To convert to approx- 
imate concentration constants subtract .09. 


b))  0.185°M 92,3=DAP.HC! yp = 0.085 — 0.400 M.A = -5.23 ppm 
8 
A3 ai =6.53) ppm. 
’ 
CyenO. eo eMOrna thine Neila) ="0 45-207 55eM A = -0.27 ppm 
273 
A3 Jao =(),22 ) ppm « 
’ 
ad) 802185 M Lysine.HCl yp = 0.4 — 0.55 M, & = -0.20 ppm 
b3 ye -0.20 ppm. 
é 


e) The linear estimate of the standard deviation of each pk or 
pK is < .03. 
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Figure 19. pH dependence of the chemical 
shifts of the, protons on the 
glycyl alpha carbon and the 
lysine epsilon carbon of gly- 
cyllysine in a 0.099 M aqueous 
solution. 
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shifts of the protons on the 
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Glycine =irnvag0. Se Meecueous 
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methods developed in Chapter III for unique resonance data 


and the constants obtained are presented in Table 14. 


C.3) Discussion 


i) oa,w-Diaminocarboxylic Acids 
The macroscopic acid dissociation constants of 

the o,w-diaminocarboxylic acids presented in Table 12 in- 
dicate that the two amino groups of all four compounds 
undergo simultaneous deprotonation. Consequently, the 
acid base chemistry of these acids is described by the 
macroscheme shown in Figure 13. 

Upon examination of the macroscopic constants in 
Table 2 and! the microscopic constants in) Tables 11) and 
13, several trends become apparent. As the number of 
intervening methylene groups increases, the acidity of 
the ammonium groups decreases, as does the separation be- 
tween the two macroscopic constants. Hay and Morris (64) 


attribute the acidity of the ammonium groups in a,u- 


diaminocarboxylic acids to the combined inductive (-1) 
effects of the nearby -NHy or -NH, and -CO, groups. As 
the w-NHG is moved further away from the a-NH3 Group, the 


SeErvong: —L effect of a mearby -NHy Group 2s reduced tor 
both ammonium groups and the observed acidity of the 


molecules decreases. 


The difference between pK, and pK, reflects the 
change in the acidity of the amino groups upon deprotona- 


tion of the other group. As can be seen from the results 
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TABLE 14 


Macroscopic and Microscopic Acid Dissociation Constants* 


of a-Glycyllysine-HCl and bysyigqly cine: Hed 


ee t®) 
Glycyllysine Lysylglycine™ 
d 
PK, 8.13 WSIS 
pK, AO sys) 1 Ob FAS) 
e 
pis Beias ye 
PK, 3 g g 
Pky 93 Oe) IM ORS 7 fal 
Phi32 a g 
a) Mixed activity-concentration constants. To convert 
LO fFaApProximate concentration constants, subtract 
BeOS 
Dime Ooo Moolycy Vivycine- HCl —=102 297 02 396 M- 
Cjpee LOS Mo lysytelycine- HCl = "0-105 — 0.210 —M- 


d) The linear estimate of the standard deviation 
was < 0202 for all’ pK and pk values Jisted. 

ee tor glycyllysine, subscrape Tedenotes the -carboxy— 
iiesecia group, 2 the glycy! ammonium group rand? 3 
the lysyl ammonium group. For lysylglycine, sub- 
script I denotes the’ carboxylic acid group, 2 the 
lysyl a-ammonium group and 3 the lysyl w-ammonium 
group. 

Blepcalculatcdeby nonlinear curveritting Of £5 q data 

to Equation 742 35 4830.00" 

g) Values obtained weré not reliable. See text for 

detaris. 
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in Tables 12 and 13, a decrease in the ratio K5/K3 reflects 


a decrease in the ratio Ky 2/123 but does not indicate a 


dSGrease ne cnesratto Kj} 5/k Since the deprotonated 


Ie 


amino group has only a moderate -I effect relative to the 


+ 
stronger “NH, GEQUD,, Ceprotonation Of the a-NH, would be 
expected tO ecduce’ the acidity of the w-NH, and vice versa. 


This reduction would be more pronounced when the two amino 
groups are in close proximity. Consequently, the differ-— 


ence between the acidities of the a-NHy in Forms & ander. 


and between the acidities of the w-NH3 in, FOrumS i sand L1 


should be greater when the two ammonium groups are in 
close proximity. The effects of the number of intervening 
bonds between the alpha and omega ammonium groups upon the 
relative acidities of the two ammonium groups in 2,3-DAP, 
ornithine and lysine are presented in Table 15. The 
attenuation of the -I effect of the -CO., , with distance, 


2 


on the w-NH, is indicated by the increasing value for 


pk from 2,3-DAP to lysine. 


ep 
In view of the previous discussion and the results 

presented in Table 15, a general consideration should be 

discussed at this time. Several authors have incorrectly 


assigned pK. to the a-ammonium group and pK, to the 


7 
w-ammonium group in the a,w-diaminocarboxylic acids (54, 
64, 65). As the macroscopic constants of simultaneously 


deprotonating functional groups are composites determined 


by the acidity of both deprotonating groups, discussion 
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TABLE 15 


Effects of the Separation of the Alpha and Omega Ammonium Groups 
on their Acidities in a,w-Diaminocarboxylic Acids 


Acid Bond pK.,-pK pk -pk pk, .-pk 
Cee pees 123612 Toe ete 
2, 3-DAP- 3 2.95 2-19 .04 
b 
ORN 5 1.76 1.59 a0 
b 
Lys 6 1.55 147 02 


a) Caleulated from constants in Tables 12 and 13. 
b)) Caleulated from constants in Tabte, 11. 
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of the relative acidities of the two groups, based on the 
Macroscopic constants, are not well founded. From EBqua- 
tions 37 and 38, it can be seen that the Je ee PK, 
describes the loss of a proton when both acidic groups are 
protonated and the larger PK, reflects the deprotonation 
of both groups after the titration of the competing func- 
tional group. Consequently, the ratio K/Ks 44 reflects 
the changes in the acidity of both groups upon deproton- 
ation of the other acidic group and does not necessarily 
indicate a wide separation in the acidity of the two 
£UnNcCtional groups. This point is clearly illustrated by 
the constants presented for EDMA in Table 4. Although 
the K,/K3 Fatio is greater than 1200, the aciditysof the 
two ammonium groups is quite Similar, as is indicated by 


the k Lato! of 2-8. Similar results are aliustrated 


We ie 
in Table 15 and have been previously obtained for 
diamines (32). 

Examination of the microscopic constants for orni— 
thine and lysine presented in Tables 11 and 13 indicates 
good agreement for the respective Pk, 5 and Pk) 53 values 
but some discrepancy in the values determined for PK} 3 
ae Pky 32° For lysine, theyCMR chemical shmiits of the 
model compounds, presented in Chapter VI, are somewhat 
anomolous and thus reduce the reliability of carbon-13 


calculations for this compound. In addition, the deter- 


mination of the microscopic acid dissociation constants 
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which describe the minor pathway (kK) 3 and Ki 32 in Figure 
Wa)\yets subject to. a lack of sensitivity due to lowstau— 
tomer concentration, as discussed in Chapter IV. As can 

be seen from the results in Table 13, the reproduciba lity 
of the values for Ki3 and K3 32 decreases as the separation 
between Kio and Kj 3 increases. The relatively low concen- 
Piatt On. or Pormmonhil (higure 7). ftoraomni thine ana lysine, 
inereases the sensitivity of the calculations to the errors 
in the chemical shift estimates and is responsible at 

least in part, for the decreased reliability of the calcu- 


lated values of pk and pk In plables.3. 
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Attempts have been made to predict the microscopic 
acid dissociation constants of diamines using model com- 
pounds. (20734)— “The effect of substitiene groups such eas 
-NHG and -NH, on the PK, Of ah, amino group are considered 
tosbe additive, thus allowing the acid dissociation, con— 
stants for’ diamines to be predicted from those form mono- 
amines and symmetric diamines. The microscopic constants 
of Symmetric diamines can be calculated directly from the 


macroscopic constants. From the following scheme, we 


can see that ky = k, and Kio os Koi: 


uly 


-/78 


“ol 


9 yy 


(lhe. et 


bab 


¥« 
“9 ee 
/ 
HFi th 
- 
i Ps 
4 


4 
iat: 


ws 


nid 


118 


+ 
HaN4CH)+ NH, 
k 
n + 
HN{CH,+,NH, HN4CH,} NH, 
k 
\ Is 
+ 
HN{CH,+ NH, 
Scheme V 


From Equations 37 and 38 we can write: 


k = Kj /2 


The difference between pK, of ethylamine and pk, of 
ethylene diamine is assumed to be due to the effect of an 
-NH, group two bonds removed while the difference between 
PK, and PK > is assumed to be due to the effect of an NH, 
group two bonds removed. By combining the pK, values of 
the appropriate monoaminocarboxylic acids and the substi- 
fuente etfLects) for amino groups: “n" bonds removed, the 
MEMeGOSCODLG acid dissociation constants Of (hes fourad, u— 


diaminocarboxyilic acids Cary be predieted., Table sb con- 


tains literature values for the acid dissociation constants 
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of four symmetric diamines (66) along with the substi- 
tuent parameters calculated from the appropriate alkyl 
amines. Table 17 contains similar data for the alpha and 
omega - monoaminocarboxylic acids. By assuming the 
substituent effects to be additive, the microscopic acid 
dissociation constants for the four a,w-diaminocarboxylic 
acids under discussion were calculated. As an example, 
PK1 5 and Pk 93 of ornithine can be predicted from the pK, 
Of Z—-aminovaleric acid and the substituent effects in 


f+ 
Table 16 for -NH, and -NH, groups four bonds removed. 


3 2 
PKj3 and PK) 53 can be predicted from the pK, of 5-amino- 
valeric acid in a similar manner. The microconstants 


estimated in this way for 2,3-DAP, 2,4-DAB, ornithine and 
lysine are given in Table 18. 

Comparison of the microscopic acid dissociation 
constants presented in Table 18 with the constants in 
Tables 11 and 13 shows that the predicted functional group 
acidities are slightly greater than the observed acidi- 
ties. Although the differences may arise due to ionic 
strength differences among the data sets used, a lack of 
additivity in substituent effects is probably the primary 
reason for the poor agreement. The lack of additivity 
may be due to the disubstitution at the alpha Carponen in 
Chapter VI it is noted that the alpha carbons of alpha 
amino acids exhibit anomolous chemical shift behavior with 


respect to the additivity of chemical Shitt (perameters 
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TABLE 17 
Acid Dissociation Constants@ for Monoaminocarboxylic 
Acids 
Monoaminocarboxylic Acid pk. 
3-aminopropanoic acid LOeZe 
4-aminobutanoic acid 10.48 
5-aminovaleric acid LORS s 
6-aminocaproic acid 10.84 
2-aminopropanoic acid AIRS) 
2-aminobutanoic acid 92606 
2-aminovaleric acid OS TA 
2-aminocaproic acid She Th 


a) Data taken from Reference 3. 
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TABLE 18 


Predzcted Macroscopic and Microscopic Acid Dissociation 
Constants for a,w-Diaminocarboxylic Acids 


273 —DAP 2,4-DAB ORNITHINE LYSINE 
Pky 5 6.49 (pk 8.68 eat his. 
Pk 3 67.6 8.74 es RS) alo eras} 
PkKj 53 ERY 9.90 107.30 LO 
PKj 35 8.70 9708 929 O49 
pK, 6.36 IO 8.65 Oreplez 


pK3 9.30 eels) 10.34 BOS ou 
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derived from amine and carboxylic acid data. Similar 
effects have also been noted by ovher authors (Ga,70) 5 
The hypothesis by Hine and Li (34) that the difference be- 


tween the substituent parameter predicted K., and the K 


3 3 
obtained by potentiometry represents the KS of a hydrogen 
bonded, monoprotic tautomer requires further substantiation. 
ii) Lysine Dipeptides 

a-glysyllysine and lysylglycine were chosen as 
models for the acid-base chemistry of the e-amino group 
of lysine in peptides and proteins. The macroscopic and 
microscopic acid dissociation constants calculated from 
chemical shift titration curves are presented in Table 14. 


Comparison of the pk values tor lysine “(Table lh) 


52 


and lysylglycine indicates an increase in the acidity of 


the o-NHy Of soe ph UtLesmupOn SUDSEALULCLON Omsine cars 


boxyl group by a peptide linkage. The greater acidity of 
the a-ammonium group of lysine as compared to monoalkyl- 
emines as ducwpartiy to the inductive effect oF tie CO, 
Group a(o47, 0/7 )2.e.Wwhen Ehe coy is converted to a peptide 
Iankage, its acid strengthening effect 1s even greater 
because of the increased inductive effect of the peptide 
bond and a decreased electrostatic attraction due to the 
increased separation of the alpha NH, and the ec A 
virtually identical increase in acidity is observed when 


thesk of glycyllysine is compared to the pK, OL glycine 


14 
(9469 from reference 67). The PK} 93 of lysylglycine shows 
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little change from the equivalent value in lysine and 
the Pk, 53 of a-glycyllysine shows only a slight increase 
in acidity. These results indicate that the e-NHy of 
lysine is not significantly increased by the incorporation 
Of, Che d=aminosand carboxylic acid groups into peptide 
bonds. 

Examination of the chemical shift titration data 
ins bagures 19 and 20 shows that the first ceprotonation 
is essentially complete prior to the onset of the second 
deprotonation. It can be seen from the values presented 
in Table 14 that pk 


and pk is PK, for the two 


fey wee 123 
peptides. These observations lead to the conclusion that 
only a very small fraction of the molecules in solution 
dissociate via the pathway I > III > IV (see Figure 13). 


As a result, the analysis of the f = pH data, of both 


Bc. 
peptides, using Equation 90, is extremely Sensitive, to 
the common resonance effects discussed in Chapter IV. 
Consequently, eye reliability of the estimates obtained 
for k33 and k1 32 from the chemical shitt data ofsgiycyil— 
lysine and lysylglycine is too low to warrant either the 
presentation of the values obtained for these constants 
Oritheir utilazation in further calculations.” “The only 


values of Ki2 and K193 reported are those obtained from 


the analysis of the f, coe DHedata. 
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CHAPTER VI 


CARBON-13 CHEMICAL SHIFT PARAMETERS FOR AMINES, CARBOXYLIC 


ACIDS AND AMINO ACIDS 


pe Introduction 


During the course of the PMR work, the instrument- 
ation required to perform Fourier transform carbon-13 
magnetic resonance became available. Proton decoupled 
carbon-13 nmr spectra have several advantages relative to 
pmr spectra, including singlet resonance patterns for each 
Magnetically non-equivalent carbon and a wider chemical 
shift range. However, the chemical shift effects, upon 
deprotonation are not as rapidly attenuated in CMR as they 
are in PMR and, for the compounds studied, all the ob- 
served carbons give rise to common resonances. 

It has been shown in previous chapters that accur- 
ate estimates of the macroscopic constants of simultan- 
eously deprotonating functional groups can be determined 
directly from common resonance chemical shift data. The 
agetermmination Of Microscopic acid dissociation Constants 
requires accurate estimates of either intermediate chem- 
tea snitt values or OL the net effect of the toeal depro— 
tonation of each functional group on the observed chemical 
Shite OL a Specific resonance. 

Fairhurst (45) used carbon-13 chemical shift-pH 


data to study the acid-base chemistry of cysteine and 
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and several related compounds. Although she was able to 
determine the macroconstants, the model compound approach 
failed to provide unambiguous estimates for the micro- 
Scopic constants. The shortcomings of the model compound 
approach were assumed to be due to the sensitivity of the 
Cc chemical shifts to conformational changes brought 
about by deprotonation. 

Due to the failure of the model compound approach 
to provide accurate estimates of intermediate chemical 
shift values, it was decided to determine if a linear 
parameterization scheme could be developed with which the 
intermediate chemical shifts could be predicted with suf- 
tTacient accurecy to calculate: microscopic acid dissocration 
constants from carbon-13 chemical shift data. 

Linear parameterization schemes have been developed 
CoO predict the cme chemical shafttsmor a variety ‘Of sys-= 
tems, including alkanes (68,69) aliphatic alcohols) (70), 
amines. (71,72), carboxylic Bee (41), and amino acids 
(43,56). The first carbon-13 parameterization scheme was 
developed for alkanes by Grant and Paul (68), and was 
suggested by the systematic changes observed in the chem- 
ical shifts of a series of homologous alkanes. Their data 
indicated that the chemical shift of a given carbon could 


be predicted by an equation of the form 


gt = bet Le Aetiee (98) 
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where 6 Us the chemical shittrofsearbonun, Alea oan 
additive parameter, aa is the number of atoms in the 7th 
position relative to carbon i and B is a constant whose 
value is quite close to the chemical shift of the methane 
carbon’ 

Lindeman and Adams (69) modified the parameter 
scheme of Grant and Paul to yield improved predictions 
for branched alkanes. Their best parameter set was based 


on the equation 


Sl aG= 9 B= 
om S > Dutsm + YSN, + oe (99) 
M=2 


where Bar Aas 


Of Carbon atoms which have M attached carbons and are 


Yor and A. are constants, Du is the number 
bonded to the kth carbon and aes and Me are the number 

Of eCoLvol aLoms three and lout DOonds —emoved. . > smepEme— 
sents the number of carbon atoms bonded to the kth carbon 
and Acw is used to describe the steric configuration of 
adjacent carbon atoms. Consequently, Lindeman and Adams 
reported four different parameter sets; Lhe parameter )set 
used in the prediction of a given chemical shift depending 
upon whether the kth carbon is a primary, secondary, 
tertiary or quaternary carbon. For substituted alkanes, 
parameterization schemes have been developed in which the 


carbon-13 chemical shift is predicted by replacing the 


substituent with CH3, CHo, or CH and then employing the 
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appropriate empirical parameters to account for the effect 
of the substituent. The validity of this approach for 
aqueous solutions has been questioned because of the dif- 
ferent environments in which the alkane shifts were mea- 
Surea (56). 

The linear parameterization scheme required for the 
prediction of the intermediate chemical shifts of poly- 
functional molecules must accurately predict the chemical 
Shifts of the carbons of interest in an’ aqueous environ— 
ment and at a given degree of protonation. To develop a 
parameterization scheme for amino acids, the carbon-13 
MME Spectra Of a total of forty-five amines, amino acids, 
and carboxylic acids were obtained. Several approaches 
to the development of the parameterization scheme have 
been investigated, including individual sets of parameters 
FOr primary, secondary and tertiary carbons, parameters 
derivedystricliy from amino aciducdata, vandla parameterssee 
which includes carbons from all compounds examined, re- 
gardless of the degree of substitution. 

As was noted in Chapter V, intermediate chemical 
shift values predicted by the linear parameterization 
schemes which are derived 2n this chapter are not on Surs 
ficient accuracy to be used in the calculation of micro- 
scopic acid dissociation constants from carbon-13 chemical 
shift data. However, the results presented in this chap- 


ter indicate that the estimates of chemical shittis obtained 
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from these parameterization schemes are more than adequate 


to aid in the assignment of carbon-13 resonances. 
DepeneouLts, and Discussion 


Dye Chemivcaie ohire Daca 
ihe Carbon—-13: chemical shifts from which the 

parameter sets are derived are presented in Tables 19-25. 
The €MR spectra of the amines and carboxylic acids were 
obtained at pH values corresponding to the fully proton- 
ated and deprotonated states. The amino acid spectra were 
measured at pH values corresponding to the protonated, 
zwitterionic, and deprotonated forms. Alkyl carbon reson- 
ance assignments were made using the effects of substi- 
tuents on the chemical shifts of directly-bonded carbon 
atoms, the dependence of the change in chemical shift 
accompanying the titration of ammonium and carboxylic acid 
Groups On distance. from the titrated group: (43,56), rela— 
tive intensities in molecules having two carbon atoms 
equivalent by symmetry and, where necessary, single fre- 
quency off-resonance decoupling. 

Carbon-13 chemical shifts have been reported prev- 
Lously for aqueous OF DO solutions of a number of the 
amines (72), carboxylic acids (41) and amino acids (40,42, 
137 oA) listed anelables 13-25 A review has also 
appeared which discusses carbon-13 studies of a wide range 


of compounds (76). Variations in experimental conditions 
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such as different temperatures, reference compounds (both 
internal and external) and concentrations make comparison 
of chemical shift data difficult. However, with few 
exceptions (73), the agreement with literature values 
obtained under conditions similar to those used in this 
thesis is within 0.5 ppm. The protonation shifts calcu— 
lated from the data in Tables 19-25 also are in agreement 
with = those reported for a variety of amines (72,75), car— 
DOs WET Crdclass( 75 )), and, aminoweac las ml 56715) 
ii) Alkyl Carbon Chemical Shift Parameters 

The carbon-13 chemical shift parameters for 
alkyl carbons were derived from the alkyl carbon chemical 
shift data in Tables 19-25 by the stepwise, multiple 


regression analysis of Equation 100. 


23 


ee Oe > aa oe (100) 
i=1 


by ts the chemical shift of the kth carbon; Sco a constant 
Which accounts for the degree of substitution, S°,;, of the 
Rea carbon; Os is an additive parameter for substituent 

ee aia N. the number Of 1, SUbStILuents:. ae is a pair 
interaction term to account for the combined effects of 
both an amino and a carboxyl group bonded to the kth car- 
bon; subscript a refers to the protonation state of the 
amino group (+ and 0 indicating protonated and deproton- 


ated respectively) and subscript c represents that of the 
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carboxyl group (0 and -). The substituent effect to which 
a parameter refers and the values derived for the param- 
eters are presented in Table 26. The parameters were 
derived from 353 observational equations, with a multiple 
Correlation coefficient cof 07998. Thesalky! carbonsor 
Malone acid was mot used! ainithe tinalsderivation of tne 
parameters because the quality of the fit for its meth- 
ylene carbon indicated the need for a pair-interaction 
parameter for carbons having two directly bonded carboxylic 
acid substituents. Insufficient data was available to 
derive such a parameter. 

The caznbon-l3° chemical shifts predicted by fhe 
parameters in Table 26 for the alkyl carbons of the mole- 
cules listed in Tables 19-25 have a standard error of the 
estimate of 0.88 ppm. Eighty-three percent of the pre- 
adietedrchemiveal shifts are: within 170" ppm of the observed 
shift and the mean of the absolute value of the residuals 
is 0.60 ppm. These results compare favourably with the 
recent work of Sarneski et al (72). These authors deter-— 
mined the carbon-13 protonation shifts of aqueous amines 
and developed a parameter set based on the alkane parameter 
set of Lindeman and Adams (69) plus five additional param- 
eters. With their parameter set, Sarneski et al predicted 
the carbon-13 chemical shifts of 48 aliphatic amines, with 
69% of the predicted shifts falling within 1 ppm of the 


observed shifts. 
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TABLE 26 


Chemical Shift Parameters for Alkyl Carbons of Amines, 


Carboxylic Acids and Amino Acids 


a) CH3 1 bond removed from the kth carbon, CH3 
moved from the kth carbon, 


b) between amino and carboxyl substituents bonded to the kth 


carbon. 


Primary carbon 


Secondary carbon 


Tertiary carbon 


CH; 
CH3 
CH; 
CH> 
CH2 
CH 
CH 
NH2 
NH2 
NH2 
NH2 
NH> 
NH,t 
NH3* 
NH37t 
NH3+ 
NH3+ 
COss 
C027 
C027 
CO2H 
CO2H 
CO2H 


NO FP WwW nN FP oO & WwW DN FP OP WDD FP DY FY DY FY WwW PL 
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pet) 


NH; '-CO2H Tae esc Eton 


NH3;t-CO2~ Interaction 


NH2-CO27— 


Interaction 


etre. 


Value 
(ppm) 
20507 
SSS 
56.56 
=A oe 
dos 

= 0.354 
=r e0 
0.00 
0.47 
See 
Y alle! 
3.00 

cot area 
Os22 
0.20 
Gezz 

oe oS 
- 2.70 
- 0.40 
0.00 
Sere 

- 2.61 
0.24 
0.00 

=" 3305 
0.00 

- 3.82 
- 5.82 
ey pee ee 


Standard 
Error 


(ppm) 


0.26 
O55 
OF | 
0.17 
0.18 
0.28 


0.39 
Omar 
0.28 
0-19 
0.20 
Onze 
0.34 
O27 
0.18 
0.20 
Om19 
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The protonation shifts given by Batchelor, Feeney 
and Roberts (75, Table II) suggest that the prediction of 
chemical shifts for amines, carboxylic acids and amino 
acids might be improved by having a separate parameter set 
LOrscach = typeron carbon) (1°) 2°)and 3°) aneach typeror 
molecule. Parameter sets have been calculated for primary, 
secondary and ere carbons from the data in Tables 
19-25 for all three types of compounds and also from the 
amino acid data only. The four parameter sets are given 
in Tables 27 and 28. Although the standard error of the 
estimate was reduced somewhat over that obtained with the 
parameters in Table 26, the improvement is small consid- 
ering the many more parameters required to fit the data. 

As a consequence of the above results, the parameters 
givens by Batchelor ict Alt (7/5) eweres re-examined -andaweve 
found to be based on the observed shifts of single com- 
pounds. However, several of the compounds chosen give rise 
to anomalous shifts. For example, the amino protonation 


Shivteswot athe CH (3.28) ppm) and the ich (0.838 ppm) =eabons 


2 
of dlycine and o-alamine, respectively are listed for the 
protonation shifts of CH. and CH carbons alpha to the amino 
group of amino acids. The values calculated from the approp- 
riate model compound data given in this chapter are 1.54 ppm 


for the alpha CH. and 1.27 ppm for the alpha methine. 


2 


The model represented by Equation 100 considers the 


deshielding of the kth alkyl carbon, due to substitution 
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Chemical Shift Parameters? 


TABLE 27 


Carbons of Amines 


-CH— 


——___ 


50709 
=4.1/ 
2218 
Oe LO 
=0 94 
1.06 
Dee ly 
=n 22 
AACA 
Pap Abe: 
=3.21 
=Oe2/ 
=e, 
127 
-1.76 
=25 9) 
-0.37 
-0.02 
1.06 
-0.64 
-0.13 
D262 
0276 
-0.80 
0.00 
1A 
ake 


In ppm vs external TMS. 
CH3 group 1 bond removed from kth carbon, CH3 group 2 bonds 
removed from kth carbon etc. 


Parameter 
Bo ' 
cop) Oe ea 
Sch. 2° 
63 CH3 3 
Cue CHo. 
rca 
S¢ CH ub 
6, CH 2 
Sg NH. 1 
59 NH2 2 
6:9 NH, 3 
511 NH, 4 
612 NH2 5 
6:3 NHS 2 
SNES yo 
615 NHS 3 
STEN ae 
S:7 NHS 5 . 
Sie) COom 9 1 
One stCOs >, 2 
SemCOsee So 
621 CO2H 1 
So2 CO2H 2 
655 CO.H 3 

ae Nee ecos a. 
4. NH3-CO27 
5 NH2-CO” 
a) 
b) 
c) 


~CH,- 


BAL DZ 
-8 .98 
Spel 
-0.36 
-4.01 
1.80 
0.00 
Ue 
9.49 
4.68 
=2 05 
0.00 
0.00 
8.48 
0.00 
-3.09 
=0. 66 
-0.24 
6.18 
-1.14 
=0715 
2.48 


—CH, 


26 


for Primary, Secondary, and Tertiary 
, Carboxylic Acids, and Amino Acids. 


.08 


0.00 
CW) 


Parameters for the interaction between amino and carboxyl 
substituents bonded to same carbon. 


Reh) 
pea 
Sol 
.14 
~D0 
89 
45 
po 
.00 
ANY 
.00 
eI 
pe HE 
43 
.00 
.60 
.00 
ee) 
736 
-47 
= ES} 
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TABLE 28 


Chemical Shift Parameters? for Primary, Secondary and 


Tertiary Alkyl Carbons of Amino Acids 


Parameter -CH- eo Sens 
Bo 62.67 51.49 19.98 
Vike oe oni STAG 0.00 
ie fel bea mers 3.60 =0463 0.00 
63 CH3; 3 0.49 0.00 2.58 
Se CH iL. -9.22 -8.48 -5.97 
eee CHoe 2 2.20 -2.43 -1.99 
d6°° CH 1 -8.25 0.00 0.00 
OAs Clo: 22 -0.71 -3.37 0.00 
So he al 0.00 ee 0.00 
Se we 3.06 | 0.66 0.00 
Sogeekiiion ane! -2.60 -2.44 -1.96 
Oya NH> i 4 -0.30 Ona 0.00 
Si2 NH2 5 -0.20 0.47 0697 
Ope) Noe 1 -1.64 0.00 0.00 
eee ae -0.30 -4.09 -4.04 
695 NHZ 3 29) ey OP -2.64 
SyeLONH, 4 ~0.35 -0.39 -0.48 
Sig (ORS al Tee. -2.55 0.00 
nice alee 0.00 -5.40 1.09 
nee COsm.1 3 -0.33 -0.32 0.00 
SoeeCODtn al 0.00 -6.04 0.00 
S21 CO2H 2 -0.31 -6.28 0.00 
foe COsH 3 -1.09 0.00 -0.31 
6, NH3CO,H™ 0.00 -3.99 0.00 
6, NH3CO2~~ 0.00 -6.13 0.00 
po NHeCOo as 0.00 -4.05 0.00 
Weeliappm wSseexternal TMS. 

b) CH; group 1 bond removed from kth carbon, CH3 group 2 
bonds removed from kth carbon, etc. 
Cjmeparameters 1or the interaction between amino and 


carboxyl substituents bonded to the same carbon. 
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at the various carbon atoms in the molecule, to be the 
sum of substituent independent and dependent parameters. 
The well known 8 effect (77), deshielding due to substit— 


ution at carbon atoms one bond removed (C, -CH, es C.-CH,X, 


9H or cone is 


represented by Oy and Sy aS a substituent independent 


where X can be CH3, CH5, CH; NH., NH,, CO 


deshielding of 6.01 ppm upon substitution of one 8 hydro- 
gen and by oF and 6 aS an additional deshielding of 


5.07 “ppl upon further substitution (to. C. -CHXY). “The 


k 
appropriate substituent parameters, for example Sg when 
x"or Yas NH,, account Or thesdependence Or tne total 
deshielding upon the nature of the substituting group. 
The shielding of the kth carbon at carbon atoms two bonds 
removea Erom the kthscarbon, the sy etftecr, (C) -C-CH, = 

C, -C-CH5X) has been factored into similar components. 
Parameters 6. and 5 indicate an increase in shielding of 
1.53 ppm upon the initial substitution and 5 and So sug- 


geceeurthner Substitution (to 'Cl-C-Cuxy) causes an 


k 
additional substituent-independent shielding of 1.52 ppm. 

In addition to the Substitution at remote carbons, 
the model represented by Equation 100 also correlates the 
chemical shift of the kth carbon to its degree of substi- 


tution and to the nature of the directly bonded substit- 


vents. The degree to which the effect of these molecular 


features is additive is indicated by the agreement between 


calculated and observed shifts. The success of this model 
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is comparable to that achieved by Lindeman and Adams for 
branched alkanes (69). 

The alpha carbons of alpha amino acids exhibit 
anomolous chemical shift behavior with respect to the 
additivity of parameters derived from amine and carboxylic 
acid data for directly-bonded amine and carboxyl substit- 
uents. This anomolous behavior was observed previously 
for the zwitterionic form of alpha amino acids (40,42) 
and the calculated alpha carbon chemical shifts, as pre- 
dicted by a model based on the alkane parameter set of 
Grant and Paul (68), were corrected by the addition of a 
6.0-ppm upfield shiitt.. The anteraction parameters of Table 
26 were derived to account for the unusual chemical, shift 
behavier of carbons bearing both an amino and aucarboxy. 
Group wes lie> Value: Of —52582 pom Lor = issan excelent 
agreement with the value of 6.0 ppm used by Sternlicht 
and co-workers (40,42) for the zwitterionic form of amino 
acids. Examination of the other interaction parameters 
Duetable: 26 indicates: that. theymagnitude of these param 
eterceds) dependent upon the protonation state Cf the 
directly bonded amino and carboxyl substituents. Litchman 
and Grant (78) also found that similar interaction param- 
eters were necessary to correlate the cmr shifts of poly- 
halomethanes. 

Further examination of the results reported in this 


chapter and in previous publications (68,69,72,79) allows 
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several conclusions to be drawn regarding the use of linear 
parameterization schemes for the prediction of carbon-13 
chemical shifts. It appears that mutual interaction be- 
tween the functional groups in polyfunctional molecules 
produces small chemical shift effects which cannot be 
predicted from model compound data, for functional groups 
which undergo simultaneous deprotonation. For example, 

the Cy carbon resonance of n-propyl amine and the Ce 
resonance of 2-aminobutyric acid and 4-aminobutyric acid 
Shite by 4.575, 63.89 “and 5.39 ppm, respectively, upon de- 
protonation of their amino groups. The presence of the 
carboxylate group obviously alters the chemical shift 
effect of an amino group two bonds removed. For simul- 
taneously deprotonating groups, the eftect of «the second 
group cannot be observed since the influence of the two 
gpoups Cannot be isolated in any pHerange= These ettects 
are particularly evident for carbons bonded to more than 
Onewrunctional group and for the Zwitterionic form of 
amino acids (40,42,79). As a consequence, parameteriza- 
tion schemes which provide excellent results when used to 
predict the chemical shifts of the model compounds from 
men the parameters were derived, may yield only rough 


estimates for the intermediate shifts of polyfunctional 


molecules. 
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iii) Carboxyl Carbon Chemical Shift Parameters 
Parameters for the prediction of the carbon-13 
chemical shifts of the carboxyl carbons of the molecules 
in Tables 21-25 were derived from the stepwise, multiple 


regression analysis of Equation 101: 


22 


°c, = — > N.6. (101) 
i=1 


“3 isa constant, with the subscript p indicating the 
protonation state of the carboxylic acid group (- if depro- 
tonated and 0 if protonated): 6 is an additive parameter 
LOLesUuDsSt i Euente., 1yand N; is the number of such substituents. 
The interaction parameter (orn used in Equation 100 is 
unnecessary as the parameters for the amino group will 
implicitly account for any interaction effects. The 
values listed in Table 29 for the parameters were derived 
from ninety-four observational equations with a multiple 
correlation coefficient of 0.995. 

Using the parameters in Table 29 to predict chem-= 
ical shifts for the carboxylic acid carbons of the mole- 
cules in Tables 21-25, the standard error of the estimate 
is 0.56 ppm and the mean of the absolute value of the 
residuals is 0.40 ppm. Ninety-eiight percent of the calcu- 
lated shifts are within 1.0 ppm of the observed shift. 

Thee eetftect’ for Carboxyl= carbons can be determined 


from aT and O. to be a 2.56 ppm downfield shift. Further 
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TABLE 29 


Chemical Shift Parameters for Carboxyl Carbons’ of 


Carboxylic Acids and Amino Acids 


Deprotonated CO, 


Protonated 


CHa 
CHa 
CH; 
ca 
Cay 
CH 
CH 


CO2H 
CO2H 
CO2H 


CH; 1 bond removed from the carboxyl Garbony CHS 


BPW DH SF WD DD OF BW DHS HD UO B&B WH NH NY FY YH S&S WW ND 


se) 


C02 


(ppm) 


u 
ih 


Value 


82.09 
The 
0.00 
0.64 
Onze 
Uae 
ZnO 
0.00 
5.08 
0.43 
2D 
2 Ou: 
D221. 
0.00 
Or29 
O26 
4.48 
25:09 
0 
0.74 
Fash 
Nba Ew 
OD) 
8290 
Zio 
1.50 


removed from the carboxyl carbon, etc. 


Standard 
Erreos 


(ppm) 


O26 


O29 
On23 
0.18 
0543 


0.58 
0.33 
0226 
Ons! 
0.35 


RSS 
0-36 
0.28 
O.26 
Oe 
O27 
0.43 
(Os PAS) 
0.26 
0.44 
0.34 
O33 


2 bonds 
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substi tution, (to C,. ~CHXY) leads to an additional de- 

shel dingMot (2652) ppm: “As was previously noted for alkyl 
carbons, the y effect is marked by an increase in shielding. 
6.5 and or indicate. an initial shift of -1.52 ppm and fur- 
ther substitution at the gamma position results in an 
additional 1.52 ppm upfield shift. These values represent 
the substituent-independent portion of the observed shift, 
with additional parameters being used to predict the 
eCfrTects of the various substituents: 

Comparison of the substituent independent fractions 
of the 8 and y effects, as predicted by the parameters in 
Paplek26 for alkyl ‘carbons and an Table 29 forn=carboxy1 
carbons, indicates that the direction of the observed 
shift is the same for both types of carbon atoms. However, 
the 8 effect is approximately two times as great for alkyl 
carbons as for carboxyl carbons. The magnitude of the 
substituent independent fraction of the y effect is approx- 
imately the same in both cases. 

Parameters 6 and 86 from table 29 andicate, that 
deprotonation of the carboxylic acid group results ina 
deshielding of the carboxyl carbon of 3.62 ppm. From 
Table 26, parameters S54 and Sig indicate a deshielding 
of 3.72 ppm for an alkyl carbon one bond removed and, from 


6 and 6 a deshielding of 1.04 ppm for an alkyl carbon 


PLE) 19” 
two bonds removed upon deprotonation of a carboxylic acid 


group in carboxylic acids or amino acids. 
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The results presented in this chapter indicate 
that linear parameterization schemes can be developed to 
predict the carbon-13 chemical shifts of alkyl and car— 
boxyl carbons of amines, carboxylic acids and amino acids 
in aqueous solution. The mutual interaction effects of 
Simultaneously deprotonation functional groups cannot be 
Predicted withesubticient accuracy fon use Insenescalcu— 
TPetion of microscopic acid dissociation constants from 
common resonance chemical shift data. The chemical shift 
parameters should, however, be of value in the assignment 
of observed chemical shifts to the carbons of peptides 
and proteins as well as amines, carboxylic acids and 


amino acids. 
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CHAPTER VII 


INTRODUCTION 


The importance of heavy metals in environmental 
pollution is well established (1). Of the heavy metals 
commonly found in the environment, lead is of qreac con- 
cern because of its high toxicity and its widespread 
occurence, both in nature (2) and in use by man (3). Due 
to thelr use as antiknock agents in fuel for internal: 
combustion engines (3), organic compounds of lead are 
especially widespread throughout the environment. 
Although much of the lead dispersed by man is eventually 
washed into natural waters (4), very little is known 
about the aqueous solution chemistry of organolead com- 
pounds. In their monograph on the organic compounds of 
lead, Shapiro and Frey (5) note that although the first 
coordination compound of organolead was reported in 1887, 
very little subsequent work has appeared in the litera- 
ture. Studies in aqueous solution are particularly 
Sparse, most reviews concentrating on reactions in non- 
aqueous solvents (6,7) or on related systems such as 
“organomercury and —tin (3). 

Recently, Wong, Chau, and Luxon (9) have reported 
the methylation of certain organic and inorganic lead 
compounds by microorganisms in lake sediments. Although 


not all of the sediment samples examined produced tetra- 
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methyllead ((CH,) ,Pb) from inorganic lead, in all cases 


the transformation of trimethyllead ((CH Pb’) acetate 


3)3 
to (CH) ,Pb was observed. Although the organometallic 
chemistry of trimethyllead has been the subject of much 
study (3,5,7,8), its aquéous solution chemistry and 
coordination chemistry have not been characterized (8). 
Such fundamental information is essential for developing 
treatments for organolead poisoning. Inhalation or 
absorption of tetraalkyl lead compounds results in the 
presence of dead in the fluids and tissues ofthe) body, 
primarily as dissolved trialkyl lead salts (10). Chela- 
tion therapy, a common een of treatment for inorganic 
ead poisoning (11), has Ehus far proven to be anestective 
in the treatment of organic lead poisoning due to the 

lack of an effective chelating agent for organic lead 
Salcs (10). 

In response to the need for information on the 
fundamental solution chemistry and coordination chemistry 
of alkyl lead salts in an aqueous environment, a programme 
directed towards characterizing the aqueous solution 
chemistry of trimethyllead and its complexes has been 
initiated at the University of Alberta. In Part II of 
this thesis, the results of a study of trimethyllead and of 
trimethyllead complexes of selected Carboxyi me vacids 1m 
an aqueous environment, is reported. The pacia-—base 


chemistry of (CH) ,Pb™ has been studied at various 
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concentration levels by both potentiometry and proton 
magnetic resonance spectroscopy. The formation constants 
of the trimethyllead complexes of five carboxylic acids 
with PK, values ranging from 3.40 to 4.65 have been deter- 
mined from pmr chemical shift data. The relationship 
between the magnitudes of the (CHa )e Pb formation constant 


and the ligand PK, is discussed. 


A. The Study of the Solution Chemistry of Metal Complexes 
by Nuclear Magnetic Resonance Spectroscopy 


The use of nuclear magnetic resonance (nmr) spec- 
troscopy to study the solution chemistry of ions and 
complexes is well established (12) and the technique has 
been applied to a variety of systems (13). The wide range 
of application of the nmr technique is due to the variety 
of ligands which contain nmr active nuclei whose chemical 
shifts are sensitive to complexation reactions. In 
addition, for organometallic ions, the chemical shift 
and, if the metal has an isotope of I - lL, the metal- 
proton coupling constant also provide information about 
the complexation reactions at the molecular level. 
.(Naturally-occurring lead contains 22.6% lead-207, whose 
I = 4). In view of the detailed information which has 
been obtained from pmr studies of the coordination 
chemistry of methylmercury (14), proton magnetic resonance 
is the principal technique being used in EnisSe study of 


trimethyllead chemistry. 


omnes ae es a sa hie ape ied 


| a aa . ae 
& y iG. eee | Py) H ota : — it io ei bus taper t ie a aa 
~ eran APL fig ae oad vada rs 
ee) ‘ ire 
LA tt a a - 
Sep toga (aA — 
ate ~ yt reat ae, Bi aT ae Se 


a ini a Tele tet ett alsa git | 263 

_ ; | &, 7 7 

, 7 . 

70 J ‘SE :WYERPS Sercngae ees des ae) oon ad (| 
a 


* 


‘ P< * 7 : Ae 
; seeder * ait “Ybasa, aad” ay DONT 
( 519 if r+. ‘s j a ay | init 2, “Vin eek * setq {MOE 


i) WUSTRYS .70, Ded | rh¥ & o2 ‘weiter a : 
ren wrt 5 Wet fi rtd hi i an Ate 
vd 


io sted) dations 1 ake snore jteddve moped. 3 


Wy Lies 
6 i Praha TOT locus! eS a4 wags bran | “ ats ie 


EA. Leteames fs), haw obit } pina hae ae TO 
=), : 1: GY ts 4° 4 an has a mit 
as - 
fyevtec” “fy See sia ore” ; Vee fe» 7 Pre os aba 0 eS ea 
wJfou2 Li esi om oily 26 a 
6 


-~ 


acchht {9S rbase, #0 *s, perkaunal tin . par Re su 


|) ape. Slaw solv aiabict bob inde we) 
: 


bi not Audit. % pre a 


we 


einga aa . 7: cs eal, “nag 


“te 


CHAPTER VIII 


EXPERIMENTAL 


De senemicals 


The carboxylic acids were of the highest grade 
commercially available and were used without further 
purification. Sodium formate was used as the source of 
formate ligand. The trimethyllead-acetic acid studies 
were carried out using trimethyllead acetate (Alfa 
Inorganics) as the source of both (CH) Pb” and acetate. 

Trimethyllead acetate was the only water soluble 
form of trimethyllead commercially available and was con- 
verted to a stock solution of trimethyllead perchlorate 
prior to use in the acid-base and complexation studies. 
An initial attempt to extract (CH) ,PbOH from benzene, 
using a method developed for triethyllead chloride (15), 
yielded only an insoluble, gelatinous mass of unknown 
composition. Lon exchange has been used to remove acetate 
from methylmercuric hydroxide solutions (16), and the 
following method was developed for the conversion of 


Z 3 
An approximately 0.25 M solution of trimethyllead 


(CH,) ,PbO CCH. to trimethyllead perchlorate. 


acetate is passed through an anion exchange column 
(Dowexsc) x) Soe dns the hydroxide form. sine eluatewise nen 
checked for acetate by examining the pmr spectrum, at 


high signal amplitude, for the acetate resonance. (In 
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basic solution, the resonance due to the protons of the 
acetate methyl group iselocated approximately 0.67 ppm 
downfield from the tert-butyl resonance of tertiary 
butanol and is detectable at concentration levels 

> 0.001 M). If a detectable concentration of acetate is 
present, the solution is passed through the column a 
second time. The ion exchange resin was regenerated with 
100 ml of 1M NaOH prior to each passage of trimethyllead 
through the column. Due to residual impurities on the 
column, the resin was changed after every 500 ml of 


trimethyllead solution. The ion exchanged (CH PbOH 


3) 3 
solution was neutralized with concentrated perchloric 
acid and was stored, tightly sealed, in the dark (15). 
Before use, the stock solution of (CH) ,PbC10, was stand- 
ardized with NaOH as described below. 

The sodium hydroxide solutions used in all phases 
of this work were prepared from a saturated solution of 
carbonate-free NaOH and doubly distilled water. Stand- 


ardization of NaOH was carried out using potassium hydro- 


gen phthalate as described by Vogel (17). 
B. pH Measurements 


All pH measurements were made at 25 + 1°C using an 
Orion Model 701 digital pH meter equipped with a standard 
glass electrode-porous ceramic junction, Saturated calo-— 


mel electrode pair. For solutions containing perchlorate 
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ion, electrical contact between the solution and the 

calomel electrode was made via a salt bridge containing 
4.0 M NaCl solution to avoid precipitation of KC10, in 
the liquid junction. The pH meter was calibrated using 


Basher “certiried standard Solutions of pH 4/00,>7200 


andiel07300:. 
Ge NMR Measurements 


The proton magnetic resonance measurements were 
obtained on a Varian A-60-D high resolution spectrometer 
eacuea probe: temperature of 25 = °C isee Chapter 111, 
Perce) rs 

Chemical shifts were measured relative to the cen- 
tral resonance of the tetramethylammonium (TMA) ion or the 
tert-butyl resonance of t-butanol, the choice of standard 
depending upon which allowed use of the narrower sweep 
width. [he chemical shiits are neported vs Pos, as des— 


cribed in Chapter II, and are the average of at least two 


scans. 
Dee Solution Preparation 


i) PMR Studies 


Solutions for the pmr studies were prepared from 
amalirquotl of stock (CH3) ,PbC10, solution and the approp-— 
riate amount of sodium perchlorate and/or carboxylic 


acid. All solutions were prepared using doubly distilled 


water. 
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For the acid-base studies, the initial sample was 
taken at.a pHeno higher than 5. Samples were taken 
approximately every 0.3 pH unit to a maximum pH between 
ian i2. 

Melee Ghalenlsiil jeisl lope qelal= (CH) Pb" -carboxylic acid 
solutions was reduced to a value about 3.0 pH units below 
the PK, OF the acid) in, Order tosobtainethe cChemicalmehiet 
Or Lhe srw ly protonated acidsand of the tree (eyeens 
ion. FOr acids with a PK, below 3.75, a Minimum po ot 
0.75 was used. Samples were taken every 0.3 of a pH unit 
to a pH value 3.0 pH units above the PK, of the acid and 
Chen every pH unit to 1a maximum pi of 11.00: 

All pH adjustments were made with concentrated 
HC1O, or carbonate-free NaOH solution. The temperature 
wasS maintained at 25 + 1°C using a water bath. 

ii) Potentiometric Studies 
The solutions used in the potentiometric studies 
were prepared from an aliquot of the stock trimethyllead 


perchlorate. Crystalline NaClO, was added to give a 


4 
Cotaleionic strength of about 0.3) Mo) ii necessary the 
initial pH was adjusted to a value below pH 3.0 with 
HC1O, prion to the beginning of the Cicrativonw 9a 55 eLeam 
of nitrogen was passed over the solution to reduce 
exposure of the solutions to atmospheric carbon dioxide. 


All titrations were carried out using carbonate-free NaOH 


and all solutions were maintained at 25 + 1° C using a 


Water bath. 
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BE. oscalculations 


The calculations were performed using a Hewlett- 


Packard Model HP-65 programmable electronic calculator. 
F. Standardization of Trimethyllead Stock Solutions 


The ion exchanged (CH) ,Pb™ stock solutions were 
standardized by potentiometric titration using a glass 
electrode. The titrations were carried out using car- 
bonate-free NaOH and the ionic strength was controlled 
using Naclo,. 

To insure complete neutralization of all the 
(CH) ,PbOH present in solution, the initial pH was re- 
duced to a value between 2 and 3 using concentrated HC1O,. 
The solution 1s thus a mixture of a strong acid (HC1O,) 
and a weak acid Richer by) and the titration yields two 
imnplectaon points. Titrationvor the added perchloric 
acid yaelds a sharp inflection point, at pH 6.0 Jand the 

(CH) 4Pb" a rather poorly defined inflection point between 
pH LO and 11. The exact position of sche (CH) Pb" end- 
point was determined using a second derivative plot of 

the titration curve (18). The titration with base con-— 
‘yerts the trimethyllead from Qe to (CH,) ,PbOH. 
Consequently, the amount of NaOH added between the first 
and second endpoints provides the concentration of 
trimethyllead in the stock solution. 


Since the possibility of dimer formation exists in 
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basic solution, (see Chapter IX), the potentiometric 
method was verified by determining the trimethyllead con- 
centration of several stock solutions using atomic absorp- 
€10n Spectroscopy . 

The atomic absorption measurements were made at 
the 283.3 nm line using a Perkin-Elmer 290-B atomic 
absorption spectrophotometer. A working curve was pre- 
pared using Pb (NO), end a 0.1 M solution of the disodium 
salt of ethylenediaminetetraacetic acid was used as a 
common matrix. Standard solutions were prepared containing 
10 to 40 Wig per ml_of lead, and. an aliquot of they tri-— 
methyllead stock solution was diluted to yield a solution 
of approximately 20-30 ug/ml of lead. Although the atomic 
absorption results tended to bewslightly higher than the 
trimethyllead concentrations obtained by potentiometry, 
the agreement is sufficient to support the potentiometric 
standardization. 

In addition, a method based upon the oxidation of 
trimethyllead to Poe) followed by an EDTA titration (19), 
has yielded results in excellent agreement with those 
obtained by pH titration. Consequently, we may safely 
assume that only (CH,) ,PbOH is present at) thevendpoint or 


the potentiometric titration. 


: r ae 


olny rake 


es: arr 
k 4 = : 
} a. oy - » a 
‘ “y Oz er : Er anaes ra (4 Ui of ai —_ 
fd ? 
a) a 


Joan Syow peewee rincod MORE 
ae : vi a 
thre ste Peer | we bs , et baw omc! . Y £ LES nmi: 
= ‘ 7 ny Ay ial Ve ; 
5 eh 20°50, pGeigrr A ‘sdtemddbctachinaven pada J qa 
» ; ae : 
oF teina 4 ‘eas, G ‘tet ~ | 0) tt 
' _ 


bP ee ke , th 4 7 
. = sl, 
— as ue + fax 
: sity tent eageimp aaa i a to, ais 
k os 5 ! my 7 a 
~ > - aT 
eet P sry’ 7 nia > re. 5 . : some a 
reiteig BiG rorya ets wavy rin, rsa 


7 

ceev ;. 7 
7, tps 

we é 


nh cage at a pat. ry fin er et oF 


ste 


$3 a al ee ee “at dips | e.- 25e" af ant, \opaae bts 


of ves i’ sand alia nol nba v 


y Vee 


¢ hl 
ARM a 


, -. H 7 
; 2 | 
ate ot 2) \ Ba Fe. eeinne =f t hserage - sh 


ba er 
j 1af“néted.. Mo) Fonda Sin pape yy bn fe a hes) wer Uae 


ix 'at paca 
> At jounomeee ie 


CHAPTER 1X 


NUCLEAR MAGNETIC RESONANCE STUDIES OF THE AQUEOUS 


SOLUTION CHEMISTRY OF TRIMETHYLLEAD 


A. The Acid-Base Chemistry of irinethyilead 


The proton magnetic resonance spectrum of trimethyl- 
lead is shown in Figure 21. The spectrum consists of a 
Singlet flanked by two less intense satellite lines. The 
central resonance is assigned to the protons of methyl 
groups bonded to isotopes of lead having a nuclear spin 
of zero while the satellites are due to methyl groups 
bonded to’ tead-207 (1 = 4) 22:63 natural abundance). ‘The 
chemical shift is given by the position of the central 
resonance; the lead-proton coupling constant by the 
separation of the satellite lines. The chemical shift of 
the methyl resonance and the lead-proton spin-spin 
coupling constant for a trimethyllead solution containing 
no coordinating ligand other than? hydroxide ion are pH 
dependent. The chemical shift of the central resonance 
Temshown as a eLunce1oneoOb pH dn, Fagune 22> fore eC. 409M 
solution of trimethyllead perchlorate (bottom curve). 
The coupling constant is shown as a function of pH in 
Figure 23. The pH dependence is due to the formation of 
trimethyllead hydroxide complex(es) as the pH is increased 


from acidic to basic. 
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Figure 21. PMR spectrum of the methyl groups of 
tramethyllead inva 0.140 M sclurion at 
approximately neutral pH. 
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Figure 227, 


pH dependence, of the methyl protons oOo. 
trimethyllead in a 0.140 M aqueous 
solution, in an aqueous solution con- 
taining 0.153 M trimethyllead and 
0.071 M pivalic acid and in an aqueous 
solution Contarning 0.102) Mitrimethy1— 
lead and (07093 "Mepivalacracida: 
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Hagure: 23. 


pH dependence of the lead 207-proton coupling 
constant in a 0.140 M aqueous solution. 
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+ _ 
(CH) 3Pb + OH 8 (CH3) ,PbOH (1) 


: 7 [ (CH) ,PbOH] 


bn ARE a Sree 87ers (2) 
[(CH,),Pb Jaoy- 


Ky is the formation constant for (CH3) ,PbOH and aoy- is 
the hydroxyl ion activity. Waters of solvation have been 
negqvuected for simplicity. 

As was shown in Chapter III, if chemical exchange 
is rapid on the nmr time scale, the observed chemical 
shift will be a weighted average of the chemical shifts 


of the various molecular forms present in solution. For 


a system described by Equation 1 we may write: 
=e DO at Pid (3) 


where 5h is the observed chemical shift of the methyl 


Ss 
resonance; P, and Hoare EHC LiEace1OnSsoL Ehestotalecon— 
centration of trimethyllead which are in the free 

((CH) Pb) and complexed ( (CH) ,PbOH) forms, respectively; 
O¢ and Se are tnerchemical Sshittceorethe frecwandecon= 
plexed: forms. if the trimethyl (lead us present only vin 
the free and 1:l Conale ede forms, me ate Py = 1 and 


Equation 3 can be rewritten as: 


SsnErTnrrErTn cen nn nnn nnn Utada 


: A 1:1 complex contains 1 molecule of trimethyllead and 
i Vagand molecule. 
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Soa anes: Pe) Stee. © (4) 
which rearranges to 
- 6 
b 
PL oF == (5) 
c iS 


To determine Ky from the chemical shift E1tracion 


data, Equation 2 is written in terms of De and Arts 


Pp 
chu 


if 1-P. KY 


(6) 


where ait is the hydrogen ion activity and hei S Ar0 ex 
Ono 

When the chemical shift titration data represented 
by the open points in Figure 22 was analysed using Equa- 
Ctons.> and 76 ,) thetlog Ky values obtained decreased 
continuously with decreasing pH. A systematic drift of 
this type indicates the presence oh a second equilibrium 
(20): 

BY analogy toythe solution equilibria ten methy |= 
mercury. (16);21,22) dimer formation is a possible second 
equilibrium. 


+ + 
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Ge [ ((CH3) Pb) 50H] re 
[(CH,) Pb" ] [ (CH,) ,PbOH] 

If a second equilibrium of the type described by 
Equation 7 is present, reducing the total concentration 
of trimethyl lead present in solution should reduce the 
concentration of dimer formed. Consequently, the values 
obtained for log K, from the chemical shift data should 
become pH independent at a sufficiently low concentration. 
As a result, pmr titrations were performed on solutions 
containing total trimethyllead concentrations of 0.069 M, 
0.034 M, 0.014 M, and 0.004 M. 

The ionic strength was maintained at about 0.3 M 
with sodium perchlorate. The chemical shift and coupling 
constant data and the log K) values calculated using the 
chemical shifts are presented in Tables 30-33 for the 
0207 Mito 0.004 -M solutionse) Due=touthe small relative 
change observed in the coupling constant values (approx- 
imately 2 Hz out of a coupling constant of 78 Hz or ~2.5%), 
the coupling constant data was not used to characterize 
tthe acid-base equilibria. 

Examination of the results in Tables 30-33 indicates 
that the drift with pH in the calculated og Ky values is 
reduced with decreasing trimethyllead concentrations. The 


essentially constant value of 4.83 obtained with the 
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TABLE. 30 


Proton Magnetic Resonance Data and Log Ky Values For 


0.069 M Trimethyllead Perchlorate 


c 
pH Chemical eoou cal a log K . 
Shift) (Hz) - 
6.72 doe 7171.00 = 
7.43 e525 20 06 
Telt, 1508 Lid $2 526 
8.13 1.488 Djina. DalA 
8.46 1.461 Fife) 508 
8282 1.420 13205 Sys enh 
OR 2 1386 78.40 4.92 
9.41 ARASSI SIS, TRIGA L 4.84 
Slo e210 78.60 4.76 
BOLO de 93 78.80 4.69 
10.46 2 78.80 4.65 
11.14 ee C10 79220 4.59 
11.40 Lye A Sy) To 00 - 
aye i= 3 Me 
b) ppm vs DSS. 
ey Uncertainty is = 1.05 uz. 


ad) Mixed activity-concentration constant, 25°C. 


TABLE) ol 


Proton Magnetic Resonance Data and Log K, Values For 


0.034 M Trimethyllead Perchlorate® 
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sel Baa os ele ee ee 
226 eros dey 0 es 
- 00 5 0 11%30 = 
eZ E524 Jie 30 = 
eek re -Le ie 3.0) 52.00 
223 L491 Hy 0 BO 
Soke) 1.464 77.80 4.98 
200 17425 73550 4.94 
229 1309 7182.30 4.86 
54 Le 858) 7 Cre3.0 4.83 
See res OM MIA AY 4.75 
<2 Ie PAB!) (ei 4.72 
alae ASS) 78.80 A792 
. 06 Lazo 79220 4.68 
04 eae) 7900 - 
47 1253 Ueks Ow - 
jee 0. SM. 
ppm vs DSS. 
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Mixed activity-concentration constant, 25°C. 
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TABLE 32 


Proton Magnetic Resonance Data and Log K, Values For 


0.014 M Trimethyllead Perchlorate® 


s (es 
ou aaa 7h Gh. log K,° 
2.90 ao 5 ALAS) _ 
4.20 a ae es) - = 
6.44 1534 Wit 0 = 
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cil) 3 60 78.25 4.87 
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b) ppm vs DSS. 

e)eoUncertainty is ~ #0505 Hz. 

d) Mixed activity-concentration constant, 25°C. 
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TABLE 33 


Proton Magnetic Resonance Data and Log K, Values For 


0.005 M Trimethyllead Perchloratea@:P 


PH Chemical shift® log ete 
7296 Pere + 
6.35 dee oul: = 
Orne Teo26 - 
dok8 D238 = 
URE Vap22 4.70 
Ses LoS 4.70 
S39 492 4.83 
B69 1.462 4.83 
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0.004 M solution is considered to be a good estimate of 
log Ki for the conditions used. 

To verify the conclusions drawn from the nmr data, 
potentiometric titrations were carried out on solutions 
with trimethyllead concentrations of 0.07 M and 0.007 M 
and ionic strength approximately 0.3 M. By neglecting 
dimer formation (Equation 7), Ky was calculated from 
potentiometric data using Equation 2 and the following 
procedure. 

For the®titration of aa solution of (CH) ,PbC10 


4 


with NaOH, the charge balance equation is: 


CHey rb | eA) ee eae ee tet kes (ele, = Ge 


CoH] can be obtained from the pH meter readings; aay 


may be neglected in -the pH range of interest. [Na‘] is 
4 ] 


from the initial concentration of added trimethyllead 


Obtained from the volume of titrant edded and [CV0 


nl when added as Naclo,, the 


contribution of the inert salt is neglected). The sodium 


cali (Since ie Salo 


and perchlorate ion concentrations are: easily corrected 

for dilution since both the initial volume and the volume 
a A + 

of titrant added are known quantities (20). [ (CH) Pb ] 


can thus be written in terms of known quantities. 


a - + 
[(CH,),Pb*] = [0H] + [(CH,),PbC10,] - [Na’} (10) 
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where [(CH,) ,PbC10,] is the total concentration of tri— 
methyllead in solution. Since only two forms of Lea mie Gi ye 


lead are being considered, we can write: 


[ (CH,) ,PbOH] Selo CHes) 


+ 
3 PbC10,] - [ (CH,) Pb ite bare 


3 
Thus all three variables on the right hand side of Equa-— 
tion 2 can be obtained from the experimental data. 

The log K) values obtained from he potentiometric 
data at the 0.07 M concentration level exhibited the same 
range of values and continuous drift with pH as was ob- 
served with the nmr data obtained at 0.069 M whereas the 
calculations using the data obtained at 0.007 M yielded 
an essentially constant log Ky value of 4.962 The reason 
for the difference between the nmr and potentiometrically 
determined log Ky values at low concentration is not 
immediately apparent but is probably due, in part, to the 
higher trimethyllead concentration used in the potentio- 
metric experiment. Also, the pH titration is more sus-— 
Geptible to carbon dioxide inteérterence) and “to errors in 
the determination of the trimethyllead concentration. 

Assuming that Equations 1 and 7 represent the 
equilibria which exist between trimethyllead and hydroxide 


jon in basic solution, the equation for the observed 


chemical shit becomes - 
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6 <0 PO P.é + P_6 (12) 


where the subscript d represents the dimeric Ort, 

((CH,) ,Pb) JOH. Examination of the trimethyllead chemical 
shift values at high pH and at various total concentra- 
tions of trimethyllead (see Tables 30-33) shows that the 
Pimiting chemical shift, at high ph, @s independent of 

the concentration level of the trimethyllead. This indic- 
ates that either os and O4 are the same, which is very 


unlikely considering the results for CH.Hg (16) or that 


3 
only (CH, ) ,PbOH is present at high pH, which is much more 
likely considering that only one equivalent of hydroxide 
is consumed in the standardization titrations (see Chapter 
ViecD) ee ASea result, 5 and on may be obtained directly 
from the chemical shift titration curves at pH's below 

6.5 and above 11.0, respectively. However, there is no 
DHeat which the dimeric species, 1s the only formpor 
trimethyllead present and, as a consequence, 64 cannot be 
measured directly. Ko and 64 have been estimated to be 


20n6 and-1l.345 ppm, respectively, by nonlincareileast 


Squares analysis of chemical shift titration data (19). 


Bo The Determination of Formation Constants of (Trimechy1— 
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lead Complexes of Selected Carboxylic Acids 


Tables 34 through 37 contain the chemical shift vs 


pH data for both the trimethyllead and ligand resonances 
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TABLE 34 


Proton Magnetic Resonance Data For The Trimethyllead- 


Pivalic Acid System* 


pH Chemical shift” J SZ 
Pivalic Acid Trimethyllead a eguuey 
1762 ia BS be} Ih Sey) UWE 
ex AOS) 1 3 eo SiS wet rah 
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b) ppm vs DSS. 
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TABLE 35 


Proton Magnetic Resonance Data for the Trimethyllead- 
Pivalic Acid System® 


pH Chemical Shire. 
Pivalie Acid Trimethyllead 
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TABLE 36 


Proton Magnetic Resonance Data for the Trimethyllead- 
Propionic Acid System? 


Chemical shift? 
[pH Trimethyllead 
er0 Ae ys} 
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TABLE 37 


Proton Magnetic Resonance Data for the Trimethyllead- 


Acetic Acid System@ 


pH Chemical shift? 
ACetVve eid Trimethyiilead 

ae, De OS) L535 
Pee e 2.005 es. 
OS) ZrO 1 aseitk 
Se 2.074 Seal 
3.46 22001 Leo 
Sei) 2.045 20035 
4.14 2.013 roel. 
4.76 Hh Sigal lees 0:0 
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Detsul 1759300 150 L 
5.68 dee 22 1 5.0/0 
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a) 402i (CH,) ,PbO,CCH, 
b) ppm vs DSS 
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ROT ely ecand aie solutions of trimethyllead and pivalic 
acid and 1:1 solutions of trimethyllead and acetic acid 
and acetyl glycine. The chemical shift vs pH data for 
the trimethyllead resonance alone for 1:1 solutions of 
trimethyllead and propionic acid and formic acid is pre- 
sented in Tables 38 and 39. 

Figure 22 presents the chemical shift of the tri— 
methyllead resonance as a function of pH for trimethyllead 
perchlorate and 1: and 2:lsobutions of trimethy llead 
and pivalic acid. Qualitatively, the effect of complexa- 
tion can be observed as the upfield displacement of the 
trimethyllead resonance in the pH region below pH 7 in 
the solutions ‘containing the carboxylic acid. VAbove pH 7, 
the hydroxyl ion competes with the carboxylate ligand for 
the trimethyllead and by a pH of 10 essentially all of the 
trimethyllead is present as a hydroxyl complex. The small 
Dlatecauein the 1d sleploty ate approximately spH) 6.0 sindicaces 
that complexation by the pivalic acid is essentially com- 
plete at this pH. This Us notean unexpected mesultesince 
the acid is almost completely deprotonated by this pH. 

The plot for trimethyllead perchlorate shows that hydroxyl 
Porolenation does not occur to an Significant degree be- 
low pH 6.5. As a consequence, it is apparent that neg- 
lecting the trimethyllead-hydroxide equilibria will have 


a 


7) 1:1 solutions contain approximately 1 mole of tri- 
methyllead per mole of ligand. 2:1 solutions: contain 
2 moles of trimethyllead per mole of ligand. 
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TABLE 38 


Proton Magnetic Resonance Data for the Trimethyllead- 
Formic Acid System@ 


Chemical shirt? 
_pH - Trimethyi lead 
1.94 093 
2.64 LE pee) 
3.00 e530 
5230 Lyo2 6 
360 sal 
Bao i 2.) 
4.3 1 253 
5.24 the aye 8) 
7.34 ee et 


a) POO Seles (CH )3 PbC10,, .0786 M sodium formate. 


b) ppm vs DSS. = 


TABLE 39 


Proton Magnetic Resonance Data for the Trimethyllead- 


Acetylglycine System@ 


pH Chemical Shift” 

ACcty i=? = akGlycinepe aie a eee 
nar # Methyl Methylene Trimethyllead 
0785 2.056 3.988 LEONE LS 
1.41 2.054 35930 eee) 
1270 2.054 Soy Nes 
2204 ZrO aS 32969 ie5S4 
he Oi 2.054 Big lsis 13 
2204 2 Oo: Say) 153.0) 
B00 2.048 Chee Ns) eee! 
S52 2.041 3.848 525 
3200 22089 3.814 a heea yeaa 
3.94 2036 360/83 ies k 
4.47 20 She ays) Dees) 
5.47 2.033 3.740 1.514 
6.60 Ze 082 3.740 ip sii 


UPd) 220: 3.742 eee 0 ls 


a) 0s o3 iM (CH3) ,PbC10,, 0.140 M acetylglycine, 
b) ppm vs DSS, 
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no effect upon the determination of the trimethyllead- 
carboxylic acid formation constants. 


The lead-207-methyl proton coupling constant 


(J 


AIL, 1 ) was measured for the mixtures in Tables 34-37 
a 


b 
and was observed to be essentially constant at 77.0 + 
0.5 HZ over the pH range 1.0-7.0. It appears, therefore, 
that VA os ih is not sufficiently sensitive to coordin- 
ation by the carboxylate ligand to be of use in the deter- 
mination of formation constants for trimethyllead- 
carboxylate complexes. The coupling constant values ob- 
served are quite similar to the value of 77.5 Hz observed 
for 30% (CH3) ,PbF, in aqueous. solution (23). 

it has ‘been previously allustrated. (16) sor 
methylmercury-carboxylic acid systems, that the chemical 
Shaite vs. pH data for the carboxylic acid) protons =canbe 
used to determine formation constants. In contrast to 
the trimethyllead resonances, below pH 7 the observed 
chemical shifts of the carboxylic acid resonances are 
affected by the protonation=deprotonation equilibriumas 
well as by complex formation. As a consequence, the 
observed chemical shift of the ligand resonance(s) is the 
weighted average of the protonated, deprotonated and com- 
plexed limiting shifts. Accounting for the additional 
equilibrium brings about a Signtticant, 1ncrease an the 
complexity of the calculations required to determine 


formation constants from the chemical shite titration data 
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syatssatd eto ar OS lead ith > | 
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and, as a result, the formation constants for the 
trimethyllead-carboxylic acid complexes studied in this 
thesis were calculated only from the observed changes in 
the chemical shifts of the trimethyllead protons. 

Although previous studies of the aqueous solution 
chemistry of the chloro-complexes of triethyllead (24,25) 
were inconclusive with regards to the preferred coordina- 
tion state of triethyllead, a study of the trimethyllead- 
diphenylthiocarbazone system (26) indicated that a 1:1 
complex was the dominant form in aqueous solution. Also, 
the tendency for (CH,) ,Pb™ to form 1:1 complexes has been 
documented for a variety of ligands in several nonaqueous 
solvents (23,27). Consequently, only 1:1 trimethyltead— 
carboxylic acid complexes were considered for the solution 
conditions used in this work. 

The formation constant of the trimethyllead com- 


plex can be defined by Equations 13 and 14: 


Gh) + RCO, 2 (CH,),PbO,CR (13) 


[ (CH. )} -PbO- CR] 


f + - 
(ich) 22bi 1 RCOr ml 


where (CH PbO.CR represents the trimethyllead-carboxylic 


yee 
acid complex and RCO, the deprotonated, uncomplexed, 


Carboxylic acid. 
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Ky can be determined from nmr chemical shift data 
by recalling: that, for a 1:1 complex, the observed chem- 
ical shift is defined by Equation 3, with 6. referring 
Eoetne Wimicing siatt of (CH3) ,PbO5CR. The: fraction of 
the trimethyllead in the complexed form is described by 


Equation 5 and the fraction in the free form by Equation 
Tae 
= 6)/(S_ - ee) (15) 


obs 


The equation for Ke can be written in terms of P, and Dead 


£ 
P . 
he = soe (en 
P_-[RCO, J 
and eo w can be determined from: 
- = 3 + 
[RCO. }] = a([TRCOjH], P[CH3) Pb i) C72) 


where [RCO,H], and GE ges Te are the total concentra- 
tions of carboxylic acid and trimethyllead respectively. 


a is the fraction of uncomplexed, deprotonated carboxylic 


acid at a given pH. 
ee Ky/ (Ky + ayt) (18) 


Substitution of Equations 5, 15, and 17 into 


Equation 16 gives: 
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sane one e) + 
t= ree [RCO,H] , - WSBT E (CHA) 4Pb 1, a (29) 
Caeet . 


Equation 19 gives Ke in terms of the total concentration 
of ligand and trimethyllead, dpe andthe chemical sntre 
da@ea. K, and 6, are the only unknowns and can be calcu- 
lated using Equation 19 by a series of successive approx- 
imations. From Figure 22, it can be seen that the pH 
range’ 9.5 to 6.0 45 the range Se: maximum (CH,) ,Pb0,CR 
concentration. By assigning the experimentally observed 
value of oe at pi 6.0 to bur a value for Ky can be 
calculated from each chemical shift-pH data point in the 
region of the titration curve over which complexation 
occurs (as determined from chemical shift vs pH plots 

such as Figure 22). As it is unlikely that the trimethyl- 
lead is totally complexed under the solution conditions 


used, See is most likely greater than on at pHTG..0 


bs 


(recall Equation 3) and as a consequence, the assigned 
value of Oe is increased by 0.0167 ppm and a new set of 


Ke values calculated. This cycle is repeated 4 or 5 times 


to obtain several sets of bord, and average Ke values. 


(6 te the same an all “datasets and sthe Ke values are 


£ 
averaged for each 6). Fach set of constants as then 
used to predict the chemical shift of the trimethyllead 


resonance at each pH value in the data set. The goodness 
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of fit between the predicted, 6 and observed chemical 


pred’ 
shifts is determined by calculating 


2% 
R = i) 
= 14 (0 aes 6 pay 


the minimum value of R indicating the best fit. The two 
successive an values between which R will be a minimum 

are thus determined and the entire calculation is repeated 
using smaller increments (.004 ppm) between successive 
estimates of oe for the range yielding a minimum for, k. 


The So and average K, values which yield the minimum 


i 
value Of R are reported in Table 40 for the: carboxylic 
acids in Tables 34-39. The’ uncertainties an the formation 
constants are the average absolute deviations from the 


reported value of K for values calculated with the 


By 
reported bu. The uncertainty in the ce values is +0.004 
Dpms Lom ally compounds. 

There are no formation constants for trimethyllead 
complexes available for comparison with the results given 
in Tablen40r. A ast eee comparison can be made with 
the previously reported formation constants for the methyl- 
Mercury complexes of the same carboxylic acids and, to 
facilitate comparison, the log Ke values of the appropriate 
methylmercury complexes reported in Reference 16 are given 


in Table 40. The values in Table 40 show that the log Ke 


of a given methylmercury complex is approximately 2.2 
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units greater than the log Ke of the equivalent trimethyl- 
lead complex. The greatest deviation occurs for formic 
acid which yields a difference of 2.3 units between the 
two log Ke values. This correlation lends some support 

to the previous use of the methylmercury-hydroxide equil- 
ibria as models for the trimethyllead-hydroxide system. 

Of greater interest, however, is the correlation 
shown in Figure 24 between the formation constant of the 
complex and the acid dissociation constant of the Ligand. 
From the plot in Figure 24 we obtain: 


pi eet elo K, + 3.0 (20) 


Equation 20) predicts that. for carboxylic acids, with PK, 


values less than 3.0, the log K, for complexes with 


ne 
trimethyllead should be less than zero. In other words, 
no stable complexes would be expected. A nmr titration 
Of a l:l solution of tErimethyl lead: monochloroacetic 
acid (pK, = 2.75) WaS percormed and the chemicalyshnirt 
data obtained indicates that no detectable complex is 
formed. 

By analogy to the methylmercury work (16), it would 


AW ata! 


be expected that the H coupling constant would be 
sensitive to complex formation. However, the coupling- 
constant data presented in Tables 30-32 for hydroxyl com- 


lexes and the even smaller changes in J observed 
p 207 ,,-1y 
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Figure 24. 
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for the carboxylic: acid complexes Pech tiadcwcne 
lead-207 proton coupling constant is not sensitive to the 
complexation of trimethyllead by the ligands studied. 
Any changes in the hybridization of lead, due to complex- 


ation, which involve the bond angles at the lead atom 


should be reflected in changes in J (23) Pe eonLer 
2075, 15 
and Drago (23) have shown, however, that J changes 
2075, 14 
ETO 1). DunzZ for (CH PDBE Son wacer 1:0 859.0 Hz wor 


Bape eed 
3)3PbClo, in hexamethylphosphoramide and that the change 


is not solely dependent upon any apparent geometrical change 


(CH 


in the planar (cH) Pb) ion. Raman spectroscopy studies 
(28) have illustrated that (CH) 3Pb" is planar in aqueous 
solutions of perchlorate salt. It is possible, then, 
that complex formation with the carboxylic acids studied 
does not bring about any significant change in the 
geometrical) configuration or in the hybridization of the 
lead ion in trimethyllead. As a result, a system of 1:1 
trimethyllead-oxalic acid was studied to determine if the 
presence of two carboxyl groups on a single ligand mole- 
cule would bring about changes in the arrangement of 
methyl groups around the lead which in) turn would be re— 
flected’ in) the Wead-207 proton coupling ‘constant. How— 
ever, in dilute solution (-014°M (CH) Pb") no detectable 
binding occurred below pH 7.0 and in more concentrated 
solution (0.14 (CH) 3Pb'), a highly insoluble precipitate 


formed above pH 2.4. 
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